(N 



O: 

in 



> 

(N 



Expansions about the gamma for the 
distribution and quantiles of a standard estimate 

by 

Christopher S. Withers 
Apphed Mathematics Group 
Industrial Research Limited 
Lower Hutt, NEW ZEALAND 



o 

^ I Saralees Nadarajah 

School of Mathematics 



University of Manchester 
Manchester M13 9PL, UK 



Abstract: We give expansions for the distribution, density, and quantiles of an estimate, 
■ building on results of Cornish, Fisher, Hill, Davis and the authors. The estimate is assumed 

c/3 . to be non-lattice with the standard expansions for its cumulants. By expanding about a 

skew variable with matched skewness, one can drastically reduce the number of terms needed 
for a given level of accuracy. The building blocks generalize the Hermite polynomials. We 
demonstrate with expansions about the gamma. 



■ Keywords: Bell polynomials; Gamma distribution; Normal distribution. 

2- 1 Introduction and summary 

The Cornish Fisher expansions due to Cornish and Fisher (1937) and Hill and Davis (1968) 
have received applications in many areas of statistics. The Cornish Fisher expansions also 
^ ■ have applications in many applied areas, including risk measures for hedge funds, margin 

setting of index futures, structural equation models, modified sudden death tests, blind 
inversion of Wiener systems, GPS positioning accuracy estimation, steady-state simulation 
analysis, blind separation of post-nonlinear mixtures, cycle time quantile estimation, estima- 
tion of the maximum average time to flower, performance of Skart, testing and evaluation, 
load flow in systems with wind generation, Value-at-Risk portfolio optimization, quantile 
mechanics, channel capacity in communications theory, economics, financial intermediation 
and physics. Three of the most recent papers applying Cornish Fisher expansions to these 
areas are: Alfredo and Arunachalam (2011), Simonato (2011) and Zhang et al. (2011). 

The aim of this paper is to develop technical tools so that the Cornish Fisher expansions 
could have wider applications. In particular, we show how one can drastically reduce the 
number of terms needed for a given level of accuracy. The building blocks involve Bell 
polynomials and Hermite polynomials. In-built routines for these polynomials are available 
in most computer algebra packages. 

Let 9 be an unknown real parameter with a non-lattice estimate having the standard 
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asymptotic cumulant expansions 

oo 

a-rie n \ r>l, (1.1) 

i=r—l 

where aio6» = &, each ario is bounded in n, and a2ie is bounded away from zero as n increases. 
We call such an estimate a standard estimate. For example, a smooth function of a sample 
mean is a standard estimate - see Withers (1983). Formulas for the leading coefficients 
were given for parametric estimates in Withers (1982) and for non-parametric estimates in 
Withers (1983, 1988). The standardized form 

Yoie = {n/a2i0f^ {e-e) (1.2) 

has cumulants expandable as 

oo 

Kr (Fow) ~ rf/^ ^riO n-\ r > 1, (1.3) 

i=i — 1 

where 

^100 = 0, Arie = arie/al{l for {r,i) 7^ (1,0). (1.4) 

Set Pnix) = Pr(y < x) for Y = Yoig, and denote its density by Pnix). Cornish and Fisher 
(1937), and Fisher and Cornish (1960) gave expansions for Pn{x) and its inverse. These can 
be re-written in the form 

00 

Pn{x) ^ P{x) - p{x) hr{x, L), (1.5) 

r=l 

oo 

P-' {Pn{x)) ^x-Y M^, L), (1.6) 

r=l 
oo 

p-1 {P{x)) ^x + J2 9r{x, L), (1.7) 

r=l 

where P = $ and p = 4> are the distribTition and density of a standard normal variable, 
denoted N ~ M{0, 1), er{x, L) is a polynomial in both x and L = (Li, L2, . . .) for e = h, f, g, 

00 

Lr = lr/r\, Ir^Y ^r,r+j-S,d , 6 = I{r > 3), (1.8) 

3=0 

and /(•) is the indicator function. However, they gave no indication of how to truncate 
these expansions for each adjusted cumulant lr- This was remedied in Withers (1984) 
which showed: 

Theorem 1.1 With notation as above, 

00 

Pn{x) « P{X) - p{x) Y n-'l'^ hrix), (1.9) 

r=l 

00 

P-' {Pn{x)) P. X - ^n-'-/' fr{x), (1.10) 
r=l 

00 

p-1 (P(x)) ^x + Y n-^/' 9r{x) (1.11) 

r=l 
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for e = h, f , g, where er{x) is a polynomial in both x and A = {Ari}, Ayi = Arte, and again, 
P = ^, p = (j). The cumulant coefficients needed for er{x), (e = h, f,g) are as follows: 



for r = 1 

for r = 2 

for r = 3 

for r = 4 

for r = 5 

for r = 6 



An 

A22 ^43; 

A12 As3 ^54; 

A23 ^44 ^65 ; 

^13 ^34 ^55 ^76; 

^24 ^45 ^66 ^87- 



Proof: The expansions (|1.9p - p.lip are obtained by substituting Ir of (|1.8p into er{x,L) of 
([LSD-dlZD, giving 



er{x,L) « ^^eH(a;)?T- * say, 



i=0 



where 



e^(x) = ^ er-2i,i{x) = er {x, Lq) + Are say, 

0<i<r/2 



Croix) = er {x,Lo) , Are = ^ er,_2j,j(x), 

l<j<r/2 



Lq — {Lqi, Lq2t . .) , Lor — lor/rl, I 



Or 



A 



r,r—5,6i 



(1.12) 
(1.13) 

(1.14) 
(1.15) 



the leading term of □ 

To apply the expansions of Theorem 1.1, one can calculate ^{x) and its inverse using 
say NAG routines GOIEAF and GOIFAF: see http:// www. nag. co. uk/ numeric/ numeri- 
caLlibraries.asp (Column 1 of Table II of Cornish and Fisher (1937, page 318) gave the main 
quantiles of N to nine decimal places.) By (jl.lip . 6 has P(x)-quantile 



0+(Wn)'/=^X]n-'-/2 5,(x), 



(1.16) 



r=0 



where goix) = x. A drawback of the method is the increasingly large number of terms 
in each er{x) as r increases. The chief contributor is the skewness coefficient /o3 = ^32) 
followed by the bias coefficient Zqi = ^11 and the second order variance coefficient I02 = ^22- 
We now show how to remove these terms. Except for sample means, E 9 and var{9) are 
unknown, but they can be approximated by truncating their cumulant expansions. Set 



J 



Srje = ^ arion \ J > r - 1, 

i=r—l 



XT -1/2 
y JK9 = 82x0 



sue], J>0,K>1. 



(1.17) 
(1.18) 



Then, 



Kr (YjKe) ~ n'-/^ ^ A 

i=r~l 



(1.19) 
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where Aff^ = for i < J, A^^ = for i < K, and the other A^^ are given in terms of 
{Ario} of (jl.4p in Section 6. (For r > 2, A;J^^ does not depend on J.) Suppose that 0326* > 0. 
(If 0320 < 0, replace (9,9) by {-9, -9).) 

Now let w be another non-lattice standard estimate with the standard cumulant expan- 
sion: 

00 

Kr {id) ^ ariw m~\ r > 1, 

i=r—l 

where aiQw = w- Since m is arbitrary up to a multiplier, we can take m = nr for some 
constant r > 0. We also assume that a32w > 0. (If 032^, < 0, replace (w, by {—w, —w).) 
By Theorem 6.1 below, 

00 

i=r—l 

where 

aJK ^ aJK _ 7JK JJK ^ r/2-iAJK (-, r,n,^ 

^ri riO riw^ riw riw \ ^ 

By Withers and Nadarajah (2012), the expansions (|1.5p - ()1.7p . (jl.9p - (|l.lip remain true for 

Pn{x) = Pr(y < x), P{x) = Pr(X < x), (1.21) 
and px{x) = p{x) is the density of X, if 

00 

Kr{Y) - Kr{X) « n''/^ ^ Arin'', (1.22) 
i=r—l 

and X, Y are non-lattice. So, these expansions hold for 

Y = YjKe, X = YjKw, Ari = Aj^ . (1.23) 

We assume that w and the cumulant coefficients known. To apply these expansions, 

we need to be able to calculate P{x) and its inverse accurately. For X linear in gamma, 
this can be done using NAG routines GOIEFF and GOIFFF. 

Now choose r so that A'^-^ = 0, that is, 

r = (A32»M32e)'. (1.24) 

This has the effect of roughly halving the number of terms in each er{x). Table 1.1 compares 
the number of terms in er{x) for different choices of J, K. The number is written as N + M, 
where is the number of terms in er{x,LQ) and M is the number of terms in Are{x) of 
p.l4p . For example, the line eo refers to the approximation 

(n/a2ie)^^^ (9 - 9^ = {m/a2iw)^^'^ {w - w) , m = nr, 

with T of ()1.24p . The four columns after the columns headed J and K, accumulate the num- 
ber of terms needed. That is, they give the number of terms needed for error 0{n~^^^^^^'^) 
when (jl.9p - ()l.lip are truncated after r + 1 terms. The final column gives the percentage 
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savings in number of terms over the Cornish-Fisher expansion, which amounts to choosing 
u) ~ A/'(0,n-i). 

For example, to calculate a quantile of 6 to 0(n^^/^) via of (jl.lip . the Cornish- 
Fisher method requires calculating 48-t-29=77 terms, but using ( J = 3, -ftT = 4) requires only 
11+7=18 terms - a saving of 77 percent. Similarly, to calculate a p-value of 9 to 0(n~^) 
via Pn of (jl.lOp . the Cornish-Fisher method requires calculating 14+2=16 terms, but using 
J = K = 2 requires only 3+1=4 terms - a saving of 75 percent. The percentage savings 
increases with r. 

[Table 1.1 about here.] 

We draw attention to the choice J = K = 1, since then ei(x) = 0, and the gamma 
approximation has error only 0{n~^), not the usual 0{n~^^^). This is the simple approxi- 
mation 

Yiw = {n/a2ieY^'^ {o - 9 - aiwn^^^ = {m/ a2iw)^^'^ [w - w - all^„m"^) , m = nr 
with T of dLMD. 

Increasing J ox K beyond those given in the table, does not decrease the number of 
terms in er{x). However, if we want to calculate er{x) up to say r = 6, then we should 
choose J, K as given on the line for r = 6, that is, J = 3, if = 4. 

We now give er{x) up to r = 6 in terms of the generalized Hermite polynomial, 

Hr = Hrx{x) = p{xy^{-Dyp{x), D = d/dx, r > 0. (1.25) 

Hr = Hrx{x) may not be a polynomial in x, but it is a polynomial in a = (ai,a2,...), 
where 

o-r = CLrX = D^a{x), a{x) = ax{x) = — Inp(x), r > 1, (1-26) 
with generating function 

oo 

Ax{x,t) = ^a^jftVr! = ax{x + t) = -lnpx{x + t). (1.27) 

r=0 

So, oi = Hi and a{x) is convex if and only if 02 > 0. The function is simpler to compute 
than Hr, (only the first two are non-zero for the normal), but expressions in terms of a are 
longer (often much longer) than expressions in terms of H. Also Hr is easily found from its 
generating function 

00 

Bx{x,t) = Y,Hrx{x)f/r\ =px{x -t)/px{x). (1.28) 

r=0 

For e = h, f,g, Theorem 1.2 gives the much simpler form for er{x) of (|1.13p . 
Theorem 1.2 With notation as above, 

er{x) = Vr + Vre 
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for e = h, f, g, where 



In particular, 



Vr — A2i-r,ill2i-r-l, 
(r+l)/2<i<r+l 

Ari = Ari/r\. (1.29) 



Vi = i/ J > 1, 

V2=^43^f2 ifK>2, 
V3=^33^^2+^54^4 J > 2, 
V4=^44^f3+^65^5 ^ > 3, 
Vs = ^34i^2 + A55H4 + ^76^f6 J > 3, 

Ve = ^45^^3 + ^66^5 + AsrHr if K > A. 



Also, 



Vre = 0ifr<3, V4e= [4^]^ c {4^) , 

Vse = [45]o e(45) + [34] 1 e(34), 

V6e = J]{[vr]o e(7r) : vr = 5^, 46, 4^} + ^ { [7r]i e(7r) : vr = 4^, 35} + [S^] ^ e (3^) 
where the e(7r), [vrjj needed for er{x) are as follows: Firstly, there is the special case 

h{ij ■■■) = so that h (1*^2*2 • • • ) = i7ii,+2i2+--i 

by {1.31 ), where H^.r = D^H^^. The other expressions needed for er{x) are: 
Forr = 4: [4^]^ = A%/21, 

f (42) =Hj- H,Hl g (42) =H7- 2H^H^ + H^hI 
Forr = 5: [45] = ^443^4, 

/(45) = Hs- H.HsH^, 5(45) = Hs - - Hj + HiH^H^, 

[34]i = ^33^443, /(34) = Hq- HiH2H^, g{34) = Hq - i?2^4 - Hi + H1H2H3. 

Forr = 6: [5^] ^ = ^54/2!, 

/ (52) =Hg- H^Hl g (52) = i?9 - 2if4^^5 + H^Hl 

[46]o =3443 ^5, 

/(46) = Hg- H^HsH^, 5(46) = Hg - H^H^ - H^H^ + HiH^^H^, 
[A\ = A%/3\, f (43) = Hu - SHiHsHj - H2HI + mlHl 
g (43) = Hn - - 3i/4^7 + SH.H^Hj + SHIh^ 

[42]^ =^443 ^44, 

[35]i = ^33^454, /(35) = Hr- ^^1/^2^4, 9(35) = H7 - H^H^ - 

[32] 2 = :433 /2!, / (32) = i/5 - HiHl g (32) = i/5 - 2H2HS + H^hI 
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Now choose w = 1, w = Gm/m^ where Gm is a gamma variable with mean m. By 

Ai^_^ = Ar,r-i,e - (r - 1)! (^32^/2)'''' , r > 4, 
Ai^ = Ai% = A„e + driAr^r-i,e, r>2,K>2, 

Aj. f_^_i = Aj^j^-^Q = Ar^r+ie + drlArrO + dr2^r,r-l,9) T > 2, K > 3, 

where 

dri = {-r/2)A23e, ^2 = (-r/2)^230 + (~2^^)^22e- 

Sections 2 to 5 deal with the general case (|1.2ip - (jl.22p . Section 2 gives simple formulas for 
H and a of (jl.25p and (jl.26p for X a gamma variable with mean m, and so for X of p.23p . 
Examples in Section 2 include the distribution and quantiles of the sample variance and the 
Studentized mean for non-normal populations. 

Section 3 re-expresses hr{x,L) using a change of notation that does away with the 
fractional coefficients in all of the papers referred to above. Fractions are eliminated by 
giving results, not in terms of {li,l2, ■ ■ but in terms of 

[r^2^^...] = (LiVn!)(47i2!)---, (1.30) 
where Lj. = lr/r\. We shall prove the following theorem. 

Theorem 1.3 With notation as above, 

hr{x,L)= ^ GrkHk-l{x), (1.31) 
fc=r,r+2,...,3r 

where 

Crk = ^ {[vr] : i G Tirk} , 
and Tirk is the set of all partitions tt = 1*^ • • • /c**^ of k such that 

S{l)ii + --- + S{k)ik = r, 

where 

S{r) = rl{r <2) + {r- 2)I{r > 3). 

For example, 

Coo = 1, Crr = Yl [l''"''^'] • (1-32) 
0<i<r/2 

The other Crk CLf^ obtained from these using 

r—i 

Cr,r+2i = '^Cjjbr-j^i (L) , I < i < r, Li = Li+2, bri{x) = Bri{x)/i\ (1.33) 
j=0 

and Bri{x) is the ordinary Bell polynomial - see Appendix A. 
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Note that ()1.32p and (|1.33p provide convenient ways to calculate hr using MAPLE, say. 
Section 4 proves the following theorem. 

Theorem 1.4 With notation as above, fr{x, L) and Qrix, L) of U.6\} - [T7?\ ) can be expressed 
in the form 

er{x, L) = Y, { W e(^) : ^ e Sr{e)} , (1.34) 

where Sr{e) is a set of partitions it of r,r + 2, . . . ,3r. The coefficients f{iT), g{Tr) are 
polynomials in H. 

Alternatively, /(vr), ^(vr) can be written as polynomials in a, as done in Appendix D. 
/(vr) are obtained via functions Cr = Crx{x) introduced by Hill and Davis (1968). We 
express these in terms of H. 

Section 5 gives fr{x) and gr{x) of (jl.lOp - p.lip . 

Section 7 extends ()1.9p and (jl.Sp to expansions for the density and other derivatives of 

Pn{x). 

Hill and Davis (1968) gave a different motivation: the distribution of the likelihood 
ratio has an expansion of the form (|1.9p . with X chi-square, or equivalently, gamma. Given 
an expansion of the form (jl.9p . they derive (jl.lOp and (jl.lip giving fr, Qr in terms of 
(/ii, /i2, . . .). Withers and Nadarajah (2012) simplified their results using Bell polynomials. 

Another such example is when X is symmetric about zero, such as A^, or Student's t. 
In this case is an even function for r even, so that is an odd function for r odd, and 
Yn = l^oiel satisfies 

oo 

Pr (Yn <x)= P{x) - p{x) n-''/i2r(x), 

where P{x) = Pr(|A| < x) = P{x) — P{—x) has density p{x) = 2p{x). 
Appendix B gives the interesting expression for Hr in terms of a, 

Hr{x) = (ai -DY 1 = {-ly'Br (-a) , r > 0, 

where is the complete Bell polynomial, and an inverse formula for in terms of H. It 
also gives the derivatives of Hr in terms of H using the functions 

hr = (ai + DY 1 = Bria), r > 0. (1.35) 

Appendix C gives /(vr) and g^n) needed in (|1.34p in terms of H. Appendix D gives the 
same but in terms of a. Appendix E gives e(7r) for e = f, g when X is a gamma variable. 
Appendix F specializes to P{x) = <I>(x) = Pr(A^ < x), the standard normal distribution, the 
choice used by Cornish and Fisher (1937) and Fisher and Cornish (1960). We give formulas 
for some /(tt), ^'(vr) that do not hold when X is non-normal. The generating function is 
BN{x,t) = e^*~* /2 = E e(^+*^)* so that Hr = Hcr, the rth Hermite polynomial, 

Herix) = (l){x)-'^{-DY()){x) = E (x + iNY , (1.36) 
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as noted by Withers (2000). Withers and Nadarajah (2011) also extended the results of 
Cornish and Fisher (1937) to general X and gave the recurrence relation 

Hr = JiHr-i, r > 1, (1.37) 

where Ji = Hi — D; its Sections 4 and 5 gave Hr for X a standardized gamma or a Student 
random variable, using V'r+i for Hr- Its Appendix A gave er{x,L) for r < 4, e = h, f, g. 

We use the notation 

[a]j = T{a + l)/T{a - j + 1) = a{a - 1) ■ ■ ■ {a - j + 1), (1.38) 
{a)j = r(Q + j - l)/r(a - 1) = a{a + l)---{a + j -1). 

When changing variables, say from Y to X = (Y — fi)/a, it is convenient to set Px{x) = 
Pr(X < x), y = fj, + ax. Then, 

px{x) = apriy), ax{x) = -Ina + ariy), 

Ax{x,t) = -Incr + Ayiy^crt), Bx{x,t) = Byiy^at), 

Hrx{x) = a'^Hrviy), arx{x) = a^'arviv), Cr+i^x{x) = a'^Crviy)- (1-39) 



2 Expansions about the gamma, with examples 

Yw needs to have a shape parameter if we want to reduce A32 to zero. Expansions about 
X a standardized gamma or were given in Section 4 of Withers and Nadarajah (2011). 
However, it is easier to evaluate first a and then H. Theorem 2.1 gives explicit formulas for 
these for gamma random variables. 

Theorem 2.1 Let G he a gamma random variable with mean m and density y™~^e^^/r(m) 
on (0, 00). So, 

Kr{G) = (r - l)!m, Arin^ = (r - ly.m^-''/^ (2.1) 
for w = G/m, w = 1. Set 

a = m-l,y = -l/y. (2.2) 



For X = G, a, Ur of \1.26\) . Hr and the generating functions 1^1. 21 ) and U.28\) are given 
by 

aciy) =y-alny + lnr(/u), AG{y,t) - ady) = t-aln(l -p), 
«rG(y) = lir = 1) + (r - ly.ay'', r > 1, 
BG{y,t) = {l + ytre', 

Hrciy) = Ni^' = 2F0 i-r, -a,y) , (2.3) 

where pFq is the generalized hypergeometric distribution (see Section 9.14 of Gradshteyn 
and Ryzhik (1965)). Also 

YjKw = {G- ii)la, (2.4) 
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where = m/s2Kw md fi = m — suwim/ S2kY^'^ ■ By hl.39\) . for X = Yjkw, 
Hr = Hrx{x) = a^'Hraiy), ar = a^'arciy) at y = fi + ax, 

r 

Cr+i,G = r! (1 + ry) + ^ Cr+i0\ r > 1. 



i=2 



In particular, 



Cr+i,r = JJ(A;a + k - I) 



k=l 



C42 = a(18a + 7), C52 = a(72a + 23), C53 = 2a(2Q + l)(24a + 11) 
C62 = 2a(600a + 163), ces = 2a (489a^ + 600a + lOl) , 
C64 = 3a(2a + 1) (lOOa^ + 113a + 32) . 



Matching skewness by (II. Sp and (12. ip gives 

All = Aug, A22 = A22e, Ar,r-1 = ^r,r-l,0 " T^~''^^{r - 1)! 

for r > 3. So, A32 = if we choose m = nr with tV2 = 2/^329, in ([23]). Then, 

/or = ^r,r-l = ^r,r-l,y - (r - 1)! (yl32y /2)''"^ 

for r > 4. In particular, A43 = ^43y — 3A'^2y/'^^ ^54 ~ ^54F ~ ^^ggy, and ^65 
15A|2y/2- 



Example 2.1 T/ie F variable is defined by -Fni,n2 = (Xni/^i)/(Xn2/^2)) where the chi- 
square variable Xni independent of Xn2- Wishart (1947) gave expansions in powers of 
n^^ , n^^ for the cumulants of Z = hiF„j_„2. Setting n = ni + n2 say, and fi = n/ui, 
it follows that 9 = Z satisfies holds with oioe = = and the non-zero Urie given by 

arr9 = [f^ + i-lYfl] (r-l)!/2. 



a. 



r,2j+r—l,i 



2j+r~l 



+ (-l)7i'''"^'' {-4y-'B,{2j+r-2)\/{2j)\, (j,r)^(0,l). 



The Bj are given by Bq = -1, Bi = 1/6, B2 = 1/30, B3 = 1/42, B^ = 1/30, B5 = 5/66. 
(This is not the current notation for the Bernoulli numbers. Apart from Bq, his Bj is what 

we call \B2j\ today.) So, 0210 = (/2 + fi)/2. Now redefine n by n = 2/^^^j^n~^, the 
harmonic mean ofni, 112. (We can do this as n is arbitrary provided that it has magnitude 
min(ni, 112).^ So, now 0210 = 1- Writing Ori = arid, the coefficients needed for er{x) are 



ei 


an 


= if2 


/i)/2,a32- (/2'-/i)/2, 




022 


= {fl + 


fl) /2,a43 = /| + /i', 


63 


ai2 


- {fl - 


fl) /6, 033 = fl - fl 054 = 3 {ft - ft) , 


64 


a23 


= (/2 + 


A) /3, 044 = 3 (/I + ft) , = 12 (/I + A') , 


65 


ai3 


= 0,034 


= fl - ft. «55 - 12 {fl - A') , are = 60 (A' - A') , 


66 


024 


= 0,045 


- 4 (/I + A') , aee = 60 {f^ + fl) , 0^7 - 360 {fl + fl 



Cornish and Fisher (1937, page 319) and Fisher and Cornish (1960, page 216) denote 
'3^2i6»/^ by <y (2. They illustrated the quantile expansion U.16\) for ni = 24, n2 = 60, P{x) = 
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0.95, giving columns 1 to 4 of Table 2.1 using P = the normal distribution. (They give 
the exact value as .26534844 ■ ■ ■ ■) 

The picture is less rosy if the degrees of freedom are small. In this case the series must 
be truncated when divergence begins, giving an upper bound to the accuracy achievable, as 
illustrated by Table 2. 2 for degrees of freedom 5 and 5. 

We can write x^/n = Gm/m, where Gm is a gamma variable with mean m = 2n. So, 
switching to ji = 2fi = n/rui and 9 = ln{Gmi / mi) —ln{Gm2 / ^^2) , where Gmi is independent 
of Gm2 ; it follows that holds with aiog = 9 = 0, and the non-zero Orie given by 

arre = [ll + {-iril] (r-l)!/2, 

a.,2,+r-M = {-ly^^ [ll'^"-^ + i-lfll''-'-'^] B,{2j + r - 2)!/(2j)!, (j, r) / (0, 1). 
For example, the leading coefficients are 

ar,r-i,9 = {r- 2)! [7r' + (-1)'^7["'] ,r>2. 

So, a2W = j2 + 71 = ^ if we redefine n as n = 1/X^j=i2 m^^ = mim2/{mi + 7712), so that 
now holds with A^g = arig. ('ji is still given by n/rui.) 

[Tables 2.1 and 2.2 about here.] 

We end with two non-parametric examples. Suppose that we have a random sample 
Xi, . . . , Xn of size n from an unknown distribution F with mean /i = n{F) and finite 
central moments fj,r = firiF). Their empirical estimates are 

n n 

fi (Fn) = X = n-'Yl {Fn) = n-' ^ {X, - X)' , 

where Fn is the empirical distribution. Set Ur = Air//^2'^^- Withers (1983), for T{F) a 
smooth functional, the cumulants of 9 = T{Fn) have an expansion of the form (jl.ip with 
9 = T{F) and the leading cumulant coefficients ario given by Theorem 3.1 there. 

Example 2.2 The distribution and quantiles of the sample variance, 9 = fi2{Fn). (After 
scaling, this is equivalent to choosing the unbiased estimate = nfi2{Fn) / {n — 1).) So, 
9 = fj,2, 0,219 = fJ-i — fJ-2- Fhe leading Ori = arie of ( fj.3|) are given by Example 4-^ of Withers 
(1983) in terms of 5 = — 1: 

ei : an = 032 = ^6 - 3^4M2 + 2^2 ~ 6^37 

62 : 022 = 4^2 - 2A*4, 043 = Ms - 4^6Ai2 + 12^4^2 - 3A*4 - 24/i5Ai3 + 96^3A*2 - 6/i2, 

63 : ai2 0, = -3m6 + 21/i4//2 - 26^2 + IS^fg, 054 ^10 - 5/i8Af2 - 40^7/^3 

-10^6Ai4 + 20/i6/X2 - 30^5 + 480^5/^3 + 360^4Ai3 + 30^4 - 60^4^2 - 1560/^3^2 + 24Ai2- 

Example 2.3 The distribution and quantiles of the Studentized mean, Yqiq = n^^ 9, where 
9 = fi2{Fn))~^^'^{X — fi) = To(F„) say. (For a normal sample, (1 — n~^)^/^yoi6» ~ tn-i, 
but otherwise it is simpler to calculate Arig for Yqiq.) So, aiog = 0, 0219 = 1, and the other 
leading Ori = ari9 = Arie of lil.3\) are given in Example 1.2 of Withers (1989b): 



ei 

62 

63 



On = -1^3/2, 032 = -2i/3, 

022 = 3 + 7j/|/2, 043 = 12 - 2vi + I2vl, 

ai2 = (—25^3 + 6t^5 — \bv^V4) /16. 



(This reference also gives the leading h^, fr, gr when X = N.) 
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3 hr{x,L) in terms of sums of partitions 

This section can be skipped on a first reading. It proves Theorem 1.3. 

Proof of Theorem 1.3 Formulas for T-Lrk are easily derived from ()1.32p and ()1.33p . Here, 
we only note that each Tirk is a distinct set of partitions of k, and 

Ur 'Hrk = Uo<j<fc/2 

is the set of all the partitions of k. 

When X = N, equation (3.1) of Withers (1984) gave the formula 

hr{x,L) = ^^Lr^ ■■■Lrj Hr^+...+rj-l{x) / jl (3.1) 

summed over j > 1, ri > 1, . . ., rj > 1, S{ri) + • • • + S{rj) = r. Its proof follows from 
the Charlier differential series - see Withers and Nadarajah (2012), so it remains valid for 
general X. 

Let us rewrite (|3.ip in the form 

3r r 
hr{x, L) =''^CrkHk-l{x), Crk = ''^Crkj, (3.2) 

k=l 3=1 

where 

Crkj = ^ri ■ ■ ■ / j\ 

is summed over ri > 1, . . . , > 1, S{ri) + • • • + S{rj) = r. 

Now suppose that {ri, . . . ,rj} consists of ii Is, ^2 2s, • • • , ik ks. The number of ways 
this can arise is the multinomial coefficient, (j^.'^.j^)- So, we can rewrite Crkj using the [ ] 
notation of (|1.3Up in the form 

summed over 

^ik = j, In H \- kik = k, S{l)ii H h S{k)ik = r. 

(For vr a partition. Hill and Davis denote [tt] by Itx-) The last constraint can be written 
ii + 2i2 + J2a=3i'^ ~ 2)^a = So, these three constraints can be written 



a=3 



is H \- ik = j - h - 12, lii H \- kik = k, h + 2^2 + ^(a - 2 

So, 

Crk = [l*^2*^] Crkiii2, (3-3) 

ii+2i2<r 

where 

Crkiii2 ~ ^ ^ [3^ • • • A:'*] 
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is summed over 

k 

'^ia = {k- r)/2 = K say, 

a=3 
k 

y~^(a — 2)ia = r — ii — 2^2 = R say. 



a=3 

By (|A.2p . Cr-fciiia = bRK{L) at Lj = Li+2. 
By ([132]), 

Cn = [1], C22 = [1'] + [2], C33 = [1'] + [12], C44 = [1^] + [1^2] + [22] , 
C55 = [1^] + [1^2] + [12^] , = [1^] + [1^2] + [1^22] + [2^] , 

and we can write p.3p as 

r 

Crk = ^ ^ ^jj^r—j,K {Lj , 
j=0 

giving ()1.33p . since -B^i^ = for i? < ET, and Cj-k = if /c — r is odd or A; < r. So, 

hr{x,L) = ^ CrkHk-l{x) 

k=r,r+2,...,3r 
r r 

i=0 j=0 

The proof is complete. □ 
We can rewrite (|1.33p as 

i 

Cr,3r-2i = ^ Cjjbr-j,r-i, r > i > 0. (3.4) 
j=0 

Using Appendix A, we obtain some special cases of the results in Theorem 1.3. 
Corollary 3.1 With notation as above, 

r-l 

Cr,r+2 = ^jj^r-ji (3-5) 
3=0 

Cr,3r = brr (L) = [3^ , (3.6) 
a,3r-2 = [1^'-'] + 6r,r-l (L) , (1) = [3'-^] , (3.7) 

a,3r-4 = C22 [3^-'] + [l3'^-'4] + br,r-2 (L) , 6.,.-2 (1) = [3^-^42] + [S^^-^ (3.8) 

where Or = bri{L) = L.r+2 for r >\, and a term with a negative power of three is discarded. 
For example, 

C24 = i4 + L1L3 = [4] + [13], C726 = [3'] , 

C746 = Le + LiL, + C22U + C33L3 = [6] + [15] + [1^4] + [24] + [1^3] + [123], 
C48 = [4^] + [35] + [134] + C22 [32] , C22 [32] = [1232] + [232] , 

^5,11 = C22 [3^] + Cn [3^4] + [342] + [325] = [1^33] + [1324] + [23^] + [34^] + [3^5] . 
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So, using the [•] notation, all the numerical coefficients of the components ofhr{x,L) are 1: 
/ii(x,L) = [l] + [3]if2, 

h2{x,L) = ([1^] + [2]) Hi + ([13] + [4])F3 + [3^] H^, 

h^{x,L) = ([l3] + [12]) H2 + ([1^3] + [23] + [14] + [5]) + ([13^] + [34]) 
+ [3'] ^8, 

h^{x,L) = ([1^] + [1^2] + [22] ) ^3 + ([1^3] + [1^4] + [123] + [15] + [24] + [6]) H, 
+ ([1^32] + [134] + [23^] + [35] + [4^]) H7 + ([13^] + [3^4]) Hg + [3^] Hu, 

and so on. 

However, it is safer to calculate Hr{x,L) in MAPLE using (fOTT) . (fL32]) . (fL33]) . ^M - 
p.Sp to avoid the chance of missing a term. 

Note how each Crk sums over a distinct set of partitions of k, and how Crk sums over 
all distinct set of partitions of k. When X = N, the 2nd term in hr{x,L), Cr,r+2Hr+i{x), 
also occurs in fr{x,L) and gr{x,L). 

To convert a term to the form given by Cornish and Fisher (1937), Fisher and Cornish 
(1960) and Withers (1984), it is only necessary to substitute 

[1^2*2 ...] = [] {1,/kiy'' /ik\ = (ri2*2 )ll q/r\ 

k=l k=l 

where r = ^ ki^ and (1*^2*^ ' ' ') = '^'/nfc=i k\^^ik\ is the partition function. 



4 Expressions for fr{x,L) and gr{x,L) in terms of (i^, H) 

We begin by giving the proof of Theorem 1.4. 

Proof of Theorem 1.4 By equation (5.6) of Withers and Nadarajah (2012), 

r 

fr{x,L) = J2i-^)'~'ckbrk{h), (4.1) 

k=l 
r 

gr{x, L) = Y^{-l)^-^Dk Kk{h) (4.2) 

k=l 

for brk of (|1.33p . where hr = hr{x, L). So, (jl.3ip . (|4.ip and (|4.2p can be written in the form 
(ll.34p . where Sr{h) = UkUrk is the set of partitions vr = {iii2 ■ ■ ■ ) such that S{ii) + S{i2) + 

• • • = r, h{7r) = Hi^+i^^ 1. If tt is a partition of k, then the coefficients /(tt), g{7r) are 

polynomials in {Hi, . . . , Hk~i}- The exception is /i(l) = /(I) = g{l) = Hq = 1. □ 

Corollary 4.1 gives some particular cases of (I4.ip - ()4.2|) and /(•), g{-), h{-). 



Corollary 4.1 With notation as above, 

gi{x, L) = fi{x, L) = hi{x, L) = [1] + [3]F2, (4.3) 
fi = hi- C2 {hih^ + hl/2) + c^h\h2/2 - Cih\/A\, 
g4 = h4- D2 {hihs + hl/2) + Dshlh2/2 - D^hf/Al. 
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Here, Ck and are the function and operator introduced by Hill and Davis (1968): 

ci = 1, Cfc+i = Kk Ck = Kk - ■ ■ Ki I, 
Di = 1, Dk+i = Dk Jfc = Ji • • • Jfc, 

where Kk = kHi + D and Jk = kHi — D. Using the expressions for Ck in terms of sl given 
in Withers and Nadarajah (2012), we obtain 

C2 = Hi, C3 = ml - H2, C4 = l5Hf - IOH1H2 + H3, 

C5 = W5Hf - W5HfH2 + l5HiHs + WH^ - H4, 

C6 = M5Hf - 1260Ff + 2imlHi + 280HiH^ - 35H2H3 - 2IH1H4 + H5. 

The coefficient of H{ in Cr+i is E N^'' . That o/F["2^2 is -{r - 1)E N^^'/S. 
Other particular cases are 

h{k) = f{k) = g{k) = Hk-i, f (r) =0fori>2, (4.4) 
/(I, k + l) = Hk- HiHk-i = -Hk.i, g{li) = 0, (4.5) 

i 

f {V-'2) = K,(H) = Y^i-iyU - m^J{il), i > 1- (4.6) 

i=i 

The last formula, (14. 6p . is just the familiar formula for Ki in terms of the non-central 
moments, m = (mi, m2, . . .): 

Ki = mi, K2 = m2 — mf, K3 = — 3m2mi + 2mf , 
K4 = 771-4 — Am^mi — 3m2 + 12m2m\ — 6mf, ■ ■ ■ 

as given by equation (3.42) of Stuart and Ord (1987) up to i = 10. The formula 

i 

3=1 

and its inverse formula, rui = Bi{K), where k = {ki,K2, • • •), were given by Withers and 
Nadarajah (2009), and in equation (2) of Comtet (1974, page 160). For example, /(1^2) = 
K3{H) = H3- 3H1H2 + 2Hf. By gS]), hi{x,L) = fi{x,L) = gi{x,L) is given by 

Siif) = Siig) = {1, 3}, /(I) = Ho = l, /(3) = H2. 

For e(7r), Sr{e) needed by (ll.34j) to compute er{x,L), 2 < r < 6, e = f, g, see Appendix C. 

5 Expansions for standardized estimates 

Theorem 5.1 gives tools for calculating fr{x) and gr{x) of ()1.10p - ()l.ll|) . 

Theorem 5.1 Expand [tt], Crk, er{x,L) in the form 

00 00 00 

W = X^Wi"-"*; Crk = ^Crkin~\ er{x,L) = ^eri(x)n"\ 

i=0 i=0 i=0 
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Then, in terms of An of lll.29\) . lor of 111.15]) . and L^r = lor/fl, [vr](L) = [vr], Crk{L) = Crk, 
we have for e = h, f, g, 

Mo = W (Lo) , CrkO = Crk (Lq) , 

[r] = L, = lr/r\, 

[r]i = Ar^r+i-S, 5 = I{r> 3), 

hri — ^ ^ CrkiHk—1. 
k=r,r+2,...,3r 

To find hro, note that when h = h = 0, Cjj = for j > 0, so that by lll.33\) . Cr,r+2i = 
hri{L), giving 

hrO = hr {x,Lq) , 

where 

r 

hr{x,L) = (L) Hr+2i-l- 

i=l 

Also, 

Ale = Ase = 0, Age = /in = A12 + ^33^2/6, 

A4e = 621, Ase = 612 + 631, Age = 622 + 641- (5.1) 

For example, for i = 1,2 and r > 3, 

i+l^r,r— 1 

and 

[23^] ^ = ^422:432:433 + ^23:432/2, 
[2^]^ = ^22:423, 

[1^4] ^ = A?iA44/2! + 7111^12^44, [1'2] ^ = ^^1^23/2 + yliiAi2:422, 

[123] 1 = ^11:422^33 + ^11^23:432 + ^12:422:432, 

[1^3^] ^ = ^2^^32^33/2 + ^11^12^32/2, [1'2] , = Al{A2z/2 + A11A12A22, 
[iH] ^ = ^3,^33/3! + ^^1^112^32/2. 
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The terms needed for Arh for r = 4, 5, 6 are 
r = 4: /i2i= ^ C2kiHk-i, 

fc=2,4,6 

C22I = [1'] 1+^23, [1^]^ - AnAi2, 

C241 = [4]i + [13]i, [4]i = ^44, [13]i = ^11^33 + A12A32, 

C'261 = [3^]]^ = ^32^33- 

r = 5 : ei2 = Ai3 + A34i72 /or e = h,f, g, 

fc=3,5,7,9 

^331 = [I'll + [12] 1, [1^]^ =^42^^12, [12]i = ^11^23+^12^22, 

C351 = [5]i + [14] 1 + [1'3]^ + [23] 1, [5]i = A55, 

[14] 1 = Aii344 + ^12^43, 

[1^3]^ = ^^1:433/2 + ^11^12:432, [23] 1 =^22^33 +3423^32, 

C37I = [13"] , + [34]i, [132] ^ = ^11^32^33 + ^12^32/2, [34]i = ^32^44 + ^33^43, 

C391 = [3^], =^32:433/2. 

r = 6 : /i22 = ^ C2fc2-fffc-i, 

fc=2,4,6 

C222 = [l']2 + [2]2, [l2]2 = ^23/2 + ^11^13, [l2]2=^24, 

C242 = [4]2 + [13]2, [412=^45, [13]2 = AllA34 + yli2333 + Al3A32, 

C262 = [33]^ = ^32^34+^433/2, 

fc=4,6,840,12 

C44i = [l1,+ [l^2]^+[22]^, [l4]^=A?iAi2/3!, 

[1^2] ^ = A11A12A22 + ^2^^23/2, [2^] ^ = ^22^23, 

C461 = [6]i + [15]i + [C22i4]i + [^33^3)1 , [6]i - Age, [15]i = ^11:455 + ^12^54, 

[C22-^4]i — ^43^^221 + [4]lC220, [4] 1 = ^44, 
[C33-^3]i — ^32C33i + ^33^330, 

C220 = A2J2! +:422, C330 = ^?i/3! + A11A22, 

Cisi = C220 [3']i + C221 [32] g + [134] 1 + [42]^ + [35] 1, 

[32] p = ^32/2!, [134]l = AiiA32A4i + v4ii]433A43 + A12A32A43, 
[42] ^ = 3443^44, [35] 1 A32A55 + A33A54, 

^4,10,1 = [13'] , + [324] ^ , [133] ^ = 3^11^32^33 + ^12^32/3!, 

[324] ^ = A32A44/2! + A33A32A43, 

C4,12,l = [3^]^ =332^33/3!. 

By h5. the terms needed for A^e for e = f,g and r = 4, 5, 6 are given by i5.2\) for r 
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and 



r = 4 : 621 = A23H1 + [4]i + [3^]^ e (S^) + [13]i e(13) 
r = 5 : 631 = ^ {[vrji e(7r) : vr € 53(6)} , 
r = 6 : 622 = ^ {M2 e(7r) : vr € 52(6)} , 
641 = ^ {[7r]i 6(7r) : vr € Si{e)} , 



where Sr{e) is given in Appendix C for e = f,g. 



Using Theorem 5.1, gi2 is given by (j5.2p . and 



Asg = 512 + 531 



where 



531 = J]{Wi5(vr), neS-M}, ^3(5) = {5, 34, 3^, 23} . 



By Table 1.1, 54 has 8+3 terms, or 2+2 terms for X12 and 2+1 terms for X23. For r < 3, 
Are does not depend on e. If /qs = 0, then 

A2h - /01^12 = A2/ = A2g = A23H1/2 + ^45^^3/4!. 

For e = f,g, 622 needs [4]2 = A45, [2]2 = A24, [3^]2 above and 

[13]2 = ^llA34 + ^12^33 + ^13^32- 

For 541 needs [7r]i for vr G S^ig) = {6, 4^, 24, 2^, 35, 3^4, 3^ 23^}: /41 needs [vrji for vr G 
54(/) = S^{g) U {15, 1^4, 1^2, 134, 13^, 123, 1^3^, l^S}. 

6 Cumulant coefficients for Yj^e of (11.181) 

Theorem 6.1 gives A'^f^ of (11.190 in terms of Arte of ()1.4p . Its proof is outhned as fohows: 
Write s = S2Ke of (|1.17p as a2wn'~^{l + e), where e = Ylf=i^ xjn^^ , xj = a2j+i,e- So, 




By dAH) 



oo 



j=k 



So, 



oo 




j=0 
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where 



depends on 



—r/2\ 

dro = l,dri = (-r/2)xi, dr2 = (-r/2)x2 + ( 2 l^i' 



XI = A22e /(I < K), X2 = A23e I{2 < K), ■■■ . 
Theorem 6.1 and Corollary 6.1 are now immediate. 

Theorem 6.1 With notation as above, 

Ai"^ = for i < J, (6.1) 

i 

Al^ = Aue (^dii= ^ dl^^-jAlje, i > J + 1, 
i=j+i 

i 

A'll^ = Arie ®dri= ^ dr^i-jArje-, r > max(2,i - 1), 

j=r-l 

where ai hi = X]j=o ^i^i-i- 
Corollary 6.1 With notation as above, 

M!j+i,e = ^i,J+i,e, ^1,5+2,0 = ^i,J+i,6»c^ii + Aij^2,e 

and for r >2, 

A'^j._lf) = Aj.^r-l,dj Aj.^0 = Aj-rd + driAf^i — i^g, 
A'l^j^l Q = j4j. g + drlAy.,^Q + dr2Ar^i — l^g. 

For r = 2 it is simpler to use 

= 0for2<i<K, 

i 

= ^ d2,i-jA2je for i > K. 

j=K+l 



The last result of the corollary follows since K2{YjKe) — 1 = S2ooe/s2Ke — 1 = ^2^0 



7 Expansions for the density of Yqiq and its derivatives 

Withers and Nadarajah (2012) gave expansions for the density and other derivatives of 
the distribution. Theorem 7.1 gives expansions under the extended Cornish and Fisher 
assumption, = 0(1) for r > 1, and under the more general assumption (jl.3p . 
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Theorem 7.1 Under the extended Cornish and Fisher assumption, If = 0(1) for r > I, 
expansions for the density ofY^iQ, pn, and its derivatives follow from 111.5]) and h3.S\) : 

oo 

i-DY pn{x)=p{x)Y,n-''^ hir{x,L), i > 0, 

T- = l 

where 

3r 

hir{x, L) = ^ CrkHk+i{x), 
k=l 

where h^i has a different meaning from that in il.l2\) - [1.14\ )- Under the more general 
assumption il.3\) . the expansions take the form 



i-DY pn{x) =p[x)Y,n~'''^ hirix), i > 0, 



r=l 

where hir{x) is hr{x) with {H^} replaced by {Hk+i+i}. For example, 
where 

Arih = Oforr = 1,2, 

^3ih = + ^33^^i+3) 

^iih = {AuAi2 + ^23) Hi+2 + (Ai2^32 + ^lU33 + Ai4)Hi+4 + ^32^33-ffi+6. 

This follows from its Charlier expansion version given in equation (2.6) of Withers and 
Nadarajah (2012), so it remains valid for general X. See equation (4.6) there for the case 
X = N. For its multivariate extension, see Section 7 of Withers and Nadarajah (2011). 

Appendix A: The ordinary Bell polynomials 

The ordinary Bell polynomial, B,j{y), is defined in terms of a sequence y = {yi,y2, ■ ■ ■),hy 

00 

S{ty = Y,Brj{y) f (A.i) 

for j = 0, 1, . . ., where 

00 

S{t) = Y,yr.f. 

r=l 

In fact, for r > 0, Brj{y) is only a function of {xi, 1 < i < r — j + 1}. For example, 

Brjiy) = for r < j, 5,o(y) = I{r = 0), 
Brj{-y) = i-iyBrjiy), 

Bri{y) = yr, B„{y) = yl, Br+l,r{y) = ryi'"^?/2, 

r-l„, I I \ „,r-2„,2 



Br+2Ay) = ^Vl 2/3 + ( 2 2/2- 
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They are tabled on page 309 of Comtet (1974) for 1 < j < r < 10. Setting 

we can write Brj and brj of (I1.33P as 

brj{y) = Brj{y)/jl = Y,[V^T'---] (A.2) 

summed over all partitions of r, that is, over lii + 2^2 + • • • = r, subject to ii + i2 + • • • = j- 
Taking the coefficient of f in S{ty+'' = S{tyS{t)^ gives the recurrence formula 

Br,j+k{y) = Brjiy) Brk{y) = ^ Baj{y)Bbk{y). 

a+b=r 

For example, k = 1 gives 

r-l 

Br,j+i{y) = ^ Baj{y)yr-a, r > j + 1, 

a=j 

and j = /c = 1 gives 

r-l 

br2{y) = Br2{y/2) = Y,yayr-a/2, r > 2, 

a=l 
( b 

'^yayr-a, if r = 26 + 1, 

a=l 

b-1 

xl/2 + '^yayr-a, if r = 26. 

a=l 

The exponential Bell polynomial, Brj{x), is the coefficient of f^/rl in S{ty /j\ when yj = 
Xj/jl. So, 

Brj{x)/r\ = Brj{y)/jl = brj{y). 

They are tabled on pages 307, 308 of Comtet (1974) for 1 < j < r < 12. 
It is easy to show that if 1^ = aWy^, = aVxr then 

Brj{Y) = a^b'' Brjiy), Brj{X) = a^b''Brj{x). (A.3) 

The multinomial expansion can be written as S{lY/r\ = ^[1*^2*2 ...] summed over all 
partitions of r. 



Appendix B: Derivatives of Hr, and Hr as a function of a 



We first show that Hr has kth derivative 

k 



r-k 



E 

i=0 



i-lybk-^Hr+i, k>0, 



(B.l) 
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where 



In particular, 



iHi + Dyi = iHi+D)b,^i. 



bo 


= 1, bi 




= AlHf 


65 


= 5\Hf 



u3 ] 



We also prove (jl.35p , giving bk in terms of a. 
Rewriting ()1.37p as 



we obtain 



Hr i — HiHf. — Hj-^i, r > 0, 



Hr k = ^ OikiHr+i, k >0, 



where 



ak+i,i = {Hi + D) aki - Uk^i-i, 
where Oki = if i < 0. The second equation follows from Hyk+i = DHfk- This gives 



But by dOZD, 



akk = {-If, ak,k^i = {-If 'kHi, 
Oik,k-i = {-if^' yAbi, hi = aio- 



Hr = (ai - D) Hr-i = (ai - 1, r > 1, 

r 

= (-l)-i?,(-a) = J](-l)^-*S„(a), r > 0. (B.2) 

1=0 



In particular, 



Hi = ai, H2 = ai - 02, Hs = af - 3aia2 + 03, 

= af — 6a^a2 + Sa^ + 4aia3 — 04, (B.3) 

H5 = af — 10afa2 + Ibaio^ + lOafa^ — 160203 — 5aia4 + 05, (B.4) 
Hq = af - 15af 02 + 45alal - I5al + 20af 03 - 60010203 + 10o| 

— 150^04 + 150204 + 601O5 — og — IO02O3 + lOofos — 10ofo2 + of. (B.5) 

So, Comtet (1974)'s table gives Hj. in terms of a to r = 12. Replacing a by —a gives 

r 

1=0 
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proving (jl.35p . Br{sL) is called the rth complete exponential Bell polynomial. So, Comtet 
(1974) 's table pages 307-308 gives br in terms of a up to r = 12. The first six are 

bo = 1, bio = 0.1, 62 = ^2 + of, 63 = as + 3aia2 + af, 
bi = + 4aia3 + 3a2 + 60^02 + a\, 

^^5 = «5 + 5aia4 + 100203 + 100^03 + 15aia2 + 10afa2 + af , 
^^6 = 06 + 60105 + I502O4 + 10o| + 150^04 + 6O01O2O3 + 15o| 
+20o?O3 + 4501O2 + 150^02 + a\. 

For example, this gives the first four derivatives of Hj. in terms of 044, 043 above and 
042 = 662, an = —463. 

Expressions for o^ in terms of H: By (|B.ip . 

ar = Hi.r-1 = J2\- ) {-'^Tbr-i-iHi+i, r > 1. 

In particular, 

Oi = 02 = li\ — H2, 03 = 2Hf — 3H1H2 + -^3, 

04 = SlHf - I2HIH2 + ^HiH^ + ?,Hl - Hi, 

05 = 4\Hf - 60HfH2 + 20HfH3 + SOHiH^ - BHiH^ - WH2H3 + H^, 

06 = 5\Hf - 360/?f /?2 + UOHfH^ - SOH^Hi + 6H1H5 -Hq + 270HfHi 
-I2OH1H2H3 - 30F| + I5H2H4. 

A simpler way to express o^ in terms of {Hi} is to apply Faa di Bruno's rule, [4i] of 
Comtet (1974), to obtain the rth derivative of o(x) = f(p{x)) at f{p) = — Inp, in terms of 
P = {Pi,P2, ■ ■ where pr = D^p{x): 

r 

Qr = ^^fjBrjip), r > 1, 

i=i 

where fj = Fj{—p)~^ at Fj = (j — 1)!, p = p{x) is the jth derivative of f{p) at p. But 
Pr = a¥Hr, where b = —1, a = p. So, by ()A.3p . -Brj(p) = p' {—^Y BrjiH). This gives the 
simple inverse formula 

r 

ar = Y.i-lY-'{j-l)lBrj{ll),r>l. 
So, Comtet (1974)'s table pages 307-308 gives o^ in terms of H up to r = 12. For example, 

O6 = -Bei + Bq2 - 2^63 + 31^64 - 4!S65 + 51^66, 

where 

^61 = Hq, Bq2 = QH1H5 + I5H2H4, ^63 = l^HfHi + 6OH1H2H3 + 15Hl 
i?64 = 20HfH3 + A5H1H2, Bq^ = 15HfH2, Bqq = Hf. 
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This is the inverse formula to ()B.2p . These two relations are essentially the relations between 
the non-central moments and the cumulants: —a = Inp, p = e~°', K = InM, M = , 
where M, K are the moment and cumulant generating functions. So, — can be identified 
with the rth cumulant, and Hr with the rth moment. 

An alternative way to express br in terms of {aj} of ()1.26p . is to set 
ei =0, Br = ^ a^'^'aj if r > 2, 6r = K — ai — Cr- 

Then, 

^2 = h = 0, dr+i = (ai + D) 6r + ra2er-i, r > 1. 

This gives 

= a\ + Cr + 5r, r < 6, 
^4 = 3a2, 

^5 = 15aia2 + 100203, 

6q = 45a^a| + 60010203 + 15o2 + 180204 + 10a|. 



Appendix C: er{x,L) for e = f,g in terms of H 

(jl.34p gives fr{x,L), gr{x,L), r < 6 in terms of certain e(7r). We now give these in terms 
of H. First consider the case I3 = 0. Let Sor{e) denote the partitions in Sr{e) needed when 
/3 = 0. For each r we first give e(7r) covered by the special cases (14. 4p . (14. Sp . (14. 6p . 

r = 2: 5o2(/) = {4,2}, /(4) - i/g, /(2) = i?i. 
= 3 : Sosif) = {5, 14, 12}, /(5) = H^, /(12) ^ H2 - Hf = -H^.i of gS]), gl]), 

r = 4: 5o4(/) = {6,42,24,22,15,1^4,122}, /(6) = //s, / (l'2) of (glD, 

/(42) =i/7-i?ii?2, 

/(24) = i/s - i/?i?3, / (22) = i/3 - i^i , /(15) = Hr-, - HiHi = -H^.^ of 

f (124) = i/5 - -ff2-ff3 - 2i?ii/4 + 2i/2i/3. 

^■ = 5: 5o5(/) = {7,45,25,16,142,124,122,125,134,132}, f{7)=He, / (1^2) of gH), 
/(45) = i/g - i?ii?3i?4, /(25) =He- HfHi, /(16) - i/g - = -H5.1 of (03]), 

/ (142) ^Hs~ HiHj - 2HiH3H^ - H2HI + 3HfHl 
/(124) ^He- H1H5 - 2HiH2H^ - H^H^ + SH^Ha, 
f (122) ^fji- H1H3 - 3H^H2 + 3i/f , 
/ (125) =He- H2Hi - 2H1H5 + 2H^H4, 

f (iH) = ffe - SHiHr, - 3H2H4 + QHlHi - + 6H1H2H3 - GHfHs. 

For r = 6, because of the increasing number of terms, we adapt the notation of Comtet 
(1974), setting 

k- l'^2'^ ■■■ = kH\'Hl^ ■■■ (C.l) 
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when giving formulas for /(vr), ^(vr). For example, 

5-1^6 = BHfHe, 15 • 1^22 = IBHfH^ 3 • 4(11) = SH^Hn- 
Using this notation, when I3 = 0, fe{x,L) is given by 
r = 6 : 

S'o6(/) ^ {8, 5^ 46, 4^ 26, 24^, 2^4, 2^, 17, 145, 125, 1^6, 1^4^ 1^22, 1^5, 1% 1*2, 1^} , 
/(8) = Hj = 7, f (1*2) of / (52) i 9 - 142, /(46) " 9 - 135, 

/ (4^) = 11 - 3 • 137 - 23^ + 3 • 1^3^, 
/(26)^7-l25, 

/ (24^) = 9 - 1^7 - 2 ■ 135 - 123^ + 3 • 1^3^ 

/ (2^4) = 7 - 2 • 1^5 - 13^ - 1^23 + 3 • l'*3, 

/(2^) = 5 - 3 • 1^3 - 1^2 + 3 • l^ 

/(17) ^ 7 - 16 = -i/e-i of (113]), 

/(145) 9 - 18 - 135 - 14^ - 234 + 3 • 1^34, 

/(125) = 7 - 16 - 1^5 - 2 • 124 + 3 • 1^4, 

/(l26) i 7- 2 -16- 25 + 2 -125, 

/ (1^42) = 9 - 2 • 18 - 27 + 2 • 1^7 - 2 ■ 135 - 2 • 142 - 4 • 234 + 12 • 1^34 
-3^ + 9■1232-12•l33^ 

/ (1^24) ^ 7 - 2 • 16 - 25 + 1^5 - 4 • 124 + 6 • 1^4 - 2 • 13^ - 2 • 2^3 + 15 • 1^23 - 12 • l'^3, 
/ (1^22) = 5 - 2 ■ 14 - 23 - 1^3 - 6 • 12^ + 21 • 1^2 - 12 • l^ 
/ (1^5) = 7 - 3 • 16 - 3 • 25 + 6 • 1^5 - 34 + 6 • 124 - 6 • 1^4, 

/ (1*4) = 7 - 4 • 16 - 6 • 25 + 12 • 1^5 - 5 • 34 + 24 • 124 - 24 • 1^4 + 8 • 13^ + 6 • 2^3 
-36 ■ 1^23 + 24 ■ I'^S. 

A second special case for ^(Tr) is 

g{2k) = Hk+i - HiHk - Hk-i {H2 - Hf) , k>l. (C.2) 
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For example, g{2A) = H5- HiH^ - i?2^3 + HIH^. When h = 0, gr{x, L) is given by (fLMIl 
as follows: 

r = 2 : SM - {4, 2}, g(4) = i/3, 5(2) = 
r = 3: 5o3(.9) = {5}, g(5) = 

r = 4: 5o4(.9) = {6,42,24,22}, .9(6) = iJg, .9(24), 5 (2^) of JOl, 
g (42) = - 2iJ3i74 + HiHl. 

r = 5 : 5o5(.9) - {25,45,7}, g{l) = He, ff(25) of jCj), 
5(45) = iJs " ifaifs - i^l + 

r = 6: S'o6(.9) = {8,52,46,43,26,242,224,23}, 5(8) - iJy - 7, 5(26) of ([021), 
3(52) =9-2-45 + 142, 
g(46) = 9 - 36 - 45 + 135, 

g (43) = 11 - 3 • 38 - 3 • 47 + 3 • 137 + 3 • 325 + 6 • 342 - 9 • 1324 - 23^ + 3 • 123^, 

g (242) = 9 - 18 - 27 + 127 - 2 ■ 36 - 2 • 45 + 4 ■ 135 + 2 • 142 + 4 • 234 - 6 • I234 
+33 -4- 1232 + 3 -1332, 

g (224) = 7 - 2 • 16 - 2 • 25 + 3 • 125 - 2 • 34 + 4 • 124 - 3 • 1^4 
+3 • 132 + 2 • 223 - 7 • l223 + 3 • 1''3, 

g (2^) = 5 - 3 ■ 14 - 3 ■ 23 + 6 ■ 123 + 6 • 122 - 10 • 1^2 + 3 • 1^ 
When h = 0, 

h2{x,L) - = f2{x,L) = g2{x,L) = [2]i?i + [4]H^. 

For the case X = A^, it is known that gr{x,L) — I{r = l)Li does not depend on Li. This 
is a key step used in of Withers (1989a, 1983) to construct parametric and non-parametric 
confidence intervals of level say 0.95 + 0{n~^/'^), given r > 1. We have seen that this 
property is also true for general X for r < 6. We now show that it is true for all r. Set 

s = Pn{y) = Pr {Yoie - Ai„ <y) = Pn{y + Xm) = P{x) 

say. Setting gQ{x) = x gives 

00 00 
Y,n-'l^gr{x) = P-\S) = Ai„ + y = X^^ + p-\s) = n-'/\ + J^n'^'/^ [gr{x)]i^^, . 

r=0 r=0 

Taking the coefficient of n"*"/^ gives gr{x,L) = I{r = l)Li + [(7r(a^)]/i=o- 
We now give the extra terms needed when ^3 7^ 0. 

r = 2: 

f (32) =H5- HiHl /(13) = i/3 - = -H2.1 of 631). 

r = 3 : 

/(34) = He- H1H2H3, 

f (33) = Hs — 3H1H2H5 — H2 + mfHl, 

f{2i) = Hi-HlH2, 

f (132) =He- H1H5 - 2H1H2H3 - Hi + SH^Hl 
f (123) = Hi- 2H1H3 - Hi + 2HfH2. 
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r = 4 : 

/(35) = 7-124, 

/ (324) = 9 - 2 • 126 - 135 - 2^3 + 3 • 1^223, 

/ (3*) = 11 - 4 • 128 - 3 • 15^ - 6 • 2^5 + 18 • 1^225 - 2^3 + 10 • 12^ - 15 • 1^2^, 

/ (23^) = 7 - 2 • 124 - 12^ - 1^5 + 3 • 1^2^, 

/(134) = 7 - 16 - 124 - 132 - 2^3 + 3 • 1^23, 

/ (13^) = 9 - 18 - 3 • 126 - 3 • 135 - 3 • 2^5 + 9 • 1^25 -4-233 

+9 • 1^2^3 + 10 • 12^ - 15 • 1^2^, 

/(123) = 5 - 14 - 1^3 - 2 • 12^ + 3 • 1^2, 

/ (1232^ = 7 - 2 • 16 - 25 + 2 • 1^5 - 2 • 124 - 2 • 13^ - 5 • 2^3 

+9 • 12^ + 12 • 1^23 - 12 • 1^2^, 

/ (1^3) = 5 - 3 • 14 - 4 • 23 + 6 • 1^3 + 6 • 12^ - 6 • 1^2. 
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r = 5 : 

/(36) = 8 - 125, 

/ (34^) = 10 - 2 • 136 - 127 - 2^3^ + 3 • 1^23^, 
/ (3^5) = 10 - 2 • 127 - 145 + 3 • 1^224 - 2^4, 

/ (3^4) i 12 - 3 ■ 129 - 138 - 3 • 156 + 9 • 1^226 - 3 • 2^6 - 2^3^ + 9 • 1^235 - 3 ■ 2^35 
+10 • 12^3 - 15 ■ 1^233, 

/ (3^) = (14) - 5 • 12(11) - 10 ■ 158 - 10 • 2^8 + 30 • 1^2^8 
-15 • 2^5^ + 45 • 1^25^ - 10 • 2^35 

+100 • 12''5 - 150 • 1^235 - 2^4 + 15 ■ 12^3 + 10 ■ 2^ - 105 • 1^2^ + 105 • 1^2^ 

/(234) = 8 - 1^6 - 125 - 134 - 12=^3 + 3 • 1^23, 

/ (23^) = 10 - 1^8 - 3 • 127 - 3 • 145 + 9 • 1^25 - 3 ■ 12^5 + 9 ■ 1^224 

-3 • 2^4 - 12^3 + 10 • 1^2" - 15 • l'*2^ 

/ (2^3) i 6 - 123 - 2 ■ 1^4 - 1^22 + 3 ■ l'^2, 

/(135) = 8 - 17 - 125 - 134 - 2^4 + 3 • 1^24, 

/ (13^4) = (10) - 19 - 2 ■ 127 - 3 • 136 - 2 • 2^6 + 6 • 1^26 - 145 

+3 • 1^35 - 235 - 2^4 + 3 • 1^224 

-3 • 2^3^ + 6 • 1^232 + 10 • 12^3 - 15 • 1^223, 

/ (13") = (12) - 1(11) - 4 • 129 - 4 • 138 - 4 ■ 2^8 + 12 • 1^28 
-6- 156- 6 -2^6 + 18-1^2^6 

-3 ■ 25^ + 9 • 1^5^ - 18 • 2^35 + 36 • 1^235 + 60 • 12^5 - 90 • 1^225 - 2^4 - 4 • 2^3^ 

+55 • 12^3 - 60 ■ 1^2^3 + 10 • 2*^ - 105 ■ 1^2^ + 105 • 1^2^, 

/ (123^) = 8 - 17 - 4 • 125 + 3 • 1^5 - 2 • 134 - 2 • 2^4 + 6 ■ 1^24 

+13 -1223 -125 -15 -1422 

+6 • 1^23 - 3 • 12^3 - 2"*, 

/ (1^34) = 8 - 2 • 17 - 26 + 2 • 1^6 - 125 - 2 • 2^4 - 3 • 134 
+6 • 1^24 - 3 • 232 + 6 • 1^3^ 
+9 • 12^3 - 12 • 1^23, 

/ (1233) ^ (10) - 2 ■ 19 - 28 + 2 ■ 1^8 - 3 ■ 127 - 6 • 136 
-6 ■ 2^6 + 18 ■ 1^26 - 3 • 145 - 9 • 235 

+18 • 1^35 + 27 • 12^5 - 36 ■ 1^25 - 4 ■ 2^4 + 9 ■ 1^224 - 12 • 2^3^ 

+18 -12232 + 74 -1233 

-90 • 1^223 + 10 • 2^ - 95 • 122^^ + 90 • l''2^ 

/ (l223) = 6 - 2 • 15 - 24 + 124 - 6 • 123 + 6 • 1^3 - 2 • 2^ + 15 • I222 - 12 • 1''2, 

/ (134) i 6 - 3 • 15 - 3 • 24 + 6 • 124 - 32 + 6 • 123 - 6 • 1^3, 

/ (1332) = 8 - 3 • 17 - 3 • 26 + 6 • l26 - 35 + 4 • 125 - 6 • 1^5 
-6 -134 -7- 224+ 18 -1224 

-12 • 232 + 18 • 1232 + 66 • 1223 - 72 • 1^23 + 9 ■ 2'' - 72 • l22'' + 60 • l'^22, 
/ (1^3) = 6 - 4 • 15 - 7 • 24 + 12 • 124 - 4 • 3^ + 32 • 123 
-24 • 1^3 + 6 ■ 2^ - 36 ■ 1^2^ + 24 ■ 1^2. 
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We skip the extra twenty five terms needed for /g when ^ 0. The extra terms needed 
for gr{x,L) in (jl.34p when I3 ^ are as follows: 

r = 2 : 

g (32) = - ^H-2Hs + H^Hl 
r = 3 : 

5(34) = He, - H2Hi - Hi + HiH2H^, 

g (3^) =Hs- 3H2H(i - 3H3H5 + 3H1H2H5 + iH^Hi + QH2HI 
— — H2 ~\~ 311^11^ . 
5(23) = Hi- HiHs - Hi + HIH2. 
r = 4 : 

5(35)^ 7- 25 - 34+ 124, 

5 (3^4) = 9 - 2 • 27 - 3 • 36 + 2 • 126 - 45 + 135 + 2^5 + 6 • 234 
+2 • 3^ - 6 • 123^ - 2^3 
-3 -1224 + 3 -12223, 

5 (3'') = 11 - 4 • 29 - 4 • 38 + 4 • 128 + 6 • 2^7 - 6 • 56 + 24 • 236 - 18 • 12^6 + 3-152 
+ 12 • 245 + 12 • 325 - 36 ■ 1235 - 10 • 2^5 + 18 • I2225 - 36 • 2234 + 24 • 12^4 
-24 • 23^ + 72 • 12232 + 17 • 2''3 - 60 • 122^3 - 10 • 12^ + 15 • 1^2"^, 

5 (232) = 7 - 16 - 3 • 25 + 125 - 2 • 34 + 4 • 124 + 2 ■ 132 + 5 • 223 
-6- 1223 - 4- 12^ + 3- 1^22. 
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r = 5 : 

g(36) = 8 - 26 - 35 + 125, 

g (342) = 10 - 28 - 3 • 37 + 127 - 2 • 46 + 2 • 136 + 2 • 235 + 2 • 24^ 
+5 ■ 3^4 - 6 • 1234 - 3 ■ 13^ - 2^3^ + 3 • 1^23^ 

g (3^5) = 10 - 2 • 28 - 2 • 37 + 2 • 127 - 46 + 2^6 - 5^ + 145 + 4 • 235 - 3 • 12^5 

+2 ■ 24^ + 2 • 3^4 - 6 • 1234 - 2^4 + 3 • 1^224, 

g (3^4) = (12) - 3 ■ 2(10) - 4 ■ 39 + 3 ■ 129 - 48 + 138 + 3 ■ 2^8 

-3 • 57+ 15 • 237 - 9 • 12^7 

-3 • 6^ + 3 • 156 + 12 • 246 + 12 • 3^6 - 4 • 2^6 - 27 • 1236 + 9 • l^i^Q 

+3 • 25^ + 9 • 345 - 9 • 1245 - 9 • 13^5 - 15 • 2^35 + 9 • 1^235 + 6 • 12^5 

-9 ■ 2^4^ - 36 ■ 23^4 + 54 ■ 12^34 + 4 • 2''4 - 15 • 1^2^4 

-6 ■ 3^ + 36 • 123^ + 13 • 2^3^ - 45 • 1^223^ - 10 ■ 12^3 + 15 • 1^233, 

g (3^) " (14) - 5 • 2(12) - 5 • 3(11) + 5 • 12(11) + 10 • 2^(10) 

-10 -59 + 40 -239 - 30- 12^9 

-10 • 68 + 10 • 158 + 20 ■ 248 + 20 • 3^8 - 60 • 1238 - 20 • 2^8 + 30 • 1^228 
+30 -257 - 90 -2^37 + 60 -1237 

+30 • 26^ + 60 • 356 - 90 • 1256 - 90 • 2^46 - 180 • 23^6 + 360 • 12^36 
+45 • 2''6 - 150 • 122^6 

+15 • 452 - 60 • 3^5 - 45 • 135^ - 45 • 2^5^ + 45 • 1^25^ - 180 ■ 2345 + 180 ■ 12^45 

+360 • 123^5 + 210 • 2^35 - 450 • 1^2^35 - 150 • 12*5 + 150 • 1^235 

+60 • 2^4^ + 360 • 2^324 - 600 • 12^34 - 51 • 2^4 + 225 • l22''4 

+120 • 23'' - 600 • 12^3^ - 225 • 2^2,'^ + 900 • 1^233^ + 315 • 12^3 - 525 • 1^2'*3 

+10 • 2' - 105 • 1^2^ + 105 • l'*2^ 

.g(234) ^ 8 - 17 - 2 • 26 + 1^6 - 2 • 35 + 2 • 125 - 42 + 4 • 134 + 2 • 2^4 - 3 • 1^24 

+3 ■ 23^ - 3 • 1^3^ - 4 • 12^3 + 3 ■ 1^23, 

g (23'"^) " (10) - 19 - 4 ■ 28 + 1^8 - 3 ■ 37 + 6 ■ 127 - 3 • 46 

+6- 136 + 9 •226-9- 1^26 

-3 • 5^ + 6 • 145 + 21 • 235 - 9 • 1^35 - 24 • 12^5 + 9 • 1^25 + 6 • 24^ + 6 • 3^4 - 36 • 1234 
-13 • 2'''4 + 27 • 1^224 - 6 ■ 13^ - 27 • 2^3^ + 36 • 1^23^ + 55 ■ 12^3 - 45 ■ 1^223 
+4 ■ 2^ - 25 • 1^2'' + 15 • l''2^ 

g (2^3) ^ 6 - 2 • 15 - 3 • 24 + 3 • 1^4 - 32 + 7 • 123 
-3 • 1^3 + 2 • 2-'^ - 7 • 1^2^ + 3 • l''2. 

We shall not give the expressions for the 11 extra (vr) needed for when Z3 7^ 0. 

Appendix D: Expansions of fr^Qr in terms of a 

Here, we give the coefficients /(tt), g(vr) needed in (|1.34|) for /,., g^, r < 4 or 5 in terms of 
a of (jl.26p . Again we first give these for the case /a = 0. 

For r = 2,3,4, /(7r)(H) has fewer/the same/greater number of terms than /(7r)(a) in 
19/4/5 cases, and g{TT){H) has fewer/the same/greater number of terms than g{TT){si) in 
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11/2/2 cases. On the other hand, as a function of x, at is generally much briefer than 
Hk- Here, we give the coefficients needed when = up to r = 5. We again adapt the 
notation of Comtet (1974), but this time in terms of a, not H: /c • 1*^2*^ • • • = ha^^^^'i ' ' ' • 
For example, 15 • VJl?^ = ISaiagfll- MAPLE gave /(vr) in terms of a as follows: 

r = 2 : /(4) = H3 = al - 80102 + 03, /(2) = Hi = oi. 
r = 3 : /(5) = H4 is given by (|R3| . 
/(14) = —80^02 + + SaiOa — 04, 
/(12) = -02. 

r = 4 : /(6) = is given by (jB.4p . 
/ (4^) = -15afa2 + 96alal - lOSoio^ 
—2040^0203 + 1050203 + 3301O3 + 6901O3 

— 35ofo4 + IO501O2O4 — 350304 + 21a^05 — 2I02O5 — ToiOg + 07, 
/(24) = — 70^^02 + I501O2 + 90^03 — IO0203 — 50104 + 05, 

/ (2^) = -30l02 + 03, 

/(15) = — 4oJo2 + I201O2 + 60^03 — IO02O3 — 4oi04 + 05, 

/ (1^4) = 601O2 + 30^03 — 9O2O3 — 3oi04 + O5, 

/(1'2) =03. 

For /s when Z3 = 0, see Appendix D. 

When ^3 = 0, the coefficients g^n) in terms of a needed for gr{x) are as follows: 

r = 2 : .g(4) = i?3 = 1^ - 3 • 12 + 3, g{2) = Hi = ai = I. 

r = 3 : .g(5) = H4 of (|R3|. 

r^4: .g(6) = i?5 of dEH), 

g (4^) = -9 • 1^2 + 72 • 1^2^ - 87 • 12^ 

+27 ■ l'*3 - 180 • 1^23 + 99 ■ 2^3 + 63 • 13^ + 99 ■ 124 

-33 • 34 - 33 • 1^4 + 21 • 1^5 - 21 • 25 - 7 • 16 + 7, 

5(24) = -3 • 1^2 + 9 • 12^ + 6 • 1^3 - 9 • 23 - 4 • 14 + 5, 

5(2^) = -12 + 3. 

r = 5 : g{7)^He of M. 

£,(45) = -12 • 1^2 + 144 • 1^2^ - 348 • 1^2^ + 96-2* 
+42 ■ 1^3 - 480 • l''^23 + 774 • 12^3 
+258 • 1^3^ - 270 • 23^ 

-64 ■ 1^4 + 396 • 1^24 - 204 • 2^4 - 268 ■ 134 + 34 • 4^ + 55 ■ 1^5 
-165 • 125 + 55 • 35 - 28 • 1^6 + 28 • 26 + 8 ■ 17 - 8, 
5(25) = -4 • l'^2 + 24 • 1^2^ - 12 • 2^ + 10 • 1^3 
-46 • 123 + 10 • 32 - 10 • 1^4 + 14 • 24 + 5 ■ 15 - 6. 

We skip gQ. The extra terms needed for /2, • • • , /s, 52, • • • , 55 when I3 ^ 3 are given in 
Appendix D. As functions of a, the Cr needed for fr in (I4.ip are 

C2 = oi, C3 = 2al + a2, C4 = Slaf + 7aia2 + 03, 
C5 = Alaf + A&a\a2 + llaios + la\ + 04, 

ce = 5!af + 324afa2 + 147a^a3 + 127aia2 + IQaia^ + 250203 + 05. 
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So, Cfe has the same number of terms as a function of H or of a. 
We now give the extra terms needed when /s / 0: 

r = 2 : 

/ (3^) = 14 • 12^ - 8 • 1^2 + 10 • 1^3 - 10 • 23 - 5 • 14 + 5, 
/(13) = -2- 12 + 3. 
r = 3: 

/(34) = -11 • 1^2 + 42 • 1^2^ - 15 • 2^ + 19 • 1^3 - 59 • 123 

+10 • 3^ - 15 • 1^4 + 15 • 24 + 2 • 15 - 6, 

/ (3-'') = 138 ■ l'^2^ - 374 • 1^23 + 104 • 2"* 

+26 • 1^3 - 500 • 1^23 + 810 • 12^3 + 280 • 1^3^ - 280 • 23^ 

-280 • 134 - 55 • 1^4 + 405 • 1^24 - 210 • 2^4 + 35-42 

+53 • 1^5 - 165 • 125 + 56-35 

-28-1^6 + 28-26 + 8-17-8, 

/(23) -5 • 1^2 + 3 • 2^ + 4 • 13 - 4, 

/ (13^) = -14 ■ 2^'* + 24 • 1^22 + 8-1^3 

-48 • 123 + 10 • 3^ - 10 • 1^4 + 15-24+15-6, 

/ (1^3) =2-2^ + 2-13-4. 
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r = 4 : 

/(35) = -14 • 1^2 + 96 • 1^2^ - 102 ■ 12^ 

+31 • 1^3 - 206 • 1^23 + 105 • 2^3 + 70 • 13^ - 102 • 3^ 

+104 • 124 - 34 • 1^4 - 35 • 34 + 21 • 1^5 - 21 • 25 - 7 • 16 + 7, 

/ (3^4) = 222 • 1^2^ - 1094 • 1^2^ + 912 ■ 12^ 

-1053 • 1^23 + 35 • 1*^3 + 3615 • 1^2^3 

-2490 • 123^ + 810 • 1^3^ - 1259 • 2^3 + 280 • 3^ 

-91 • 1^4 + 1185 • 1^24 - 1860 • 12^4 - 1255 • 1^34 + 1260 • 234 + 315 ■ 14^ 

+121 • 1*5 - 749 • 1^25 + 503 • 135 + 378 • 2^5 - 126 • 45 

-82 • 1^6 + 250 • 126 - 84 • 36 + 36 • 1^7 - 36 • 27 - 9 • 18 + 9, 

/ (3*) = -3792 • 1^2^ + 14512 • 1^2* - 9892 • 12^ 

-1792 • 1*^23 + 27724 • 1''223 + 3200 • 1^3^ 

-43360 • 1^23^ - 64976 • 1^233 + 15400 • 1^3^ + 67910 • 12^3^ 

+17264 • 2*3 - 15400 • 23^ 

-80 • 1^4 - 31680 • 1^224 + 4760 • 1^24 - 10130 • 1*34 

+67880 • 1^234 + 5560 • 1^4^ - 17185 • 124^ + 33780 • 12^4 - 23100 • 13^4 - 34650 • 2^ 
+5775 • 34^ + 244 • 1^5 - 5992 • 1*25 + 20028 • 1^225 - 27436 • 1235 + 8956 • 1^35 

-6900 • 1^45 - 6924 • 2^5 + 6930 • 245 + 4620 • 3^5 + 1383 • 15^ 

-350 • 1^6 + 4396 • 1^26 - 6818 • 12^6 - 4620 • 1^36 + 2310 • 146 + 4620 • 236 -462-5 

+298 • 1*7 - 1980 • 1^27 + 32 • 1^27 + 1320 • 137 + 990 • 2^7 - 330 • 47 

+491 • 128 - 161 • 1^8 - 165 • 38 + 55 • 1^9 -55-29 

-11 -1(10) + (11), 

/ (23^) = 72 • 1^2^ - 98 • 12^ 

+17 • 1*3 - 192 • 1^23 + 70 • 13^ + 105 • 2^3 - 28 • 1^4 

+103 • 124 - 35 • 34 + 20 • 1^5 - 21 • 25 - 7 • 16 + 7, 

/(134) = -84 • 12^ - 141 • 1^23 + 44 • 1^2^ + 11 • 1*3 + 59 • 13^ + 104 • 2^3 - 19 • 1^4 

+89 • 124 - 35 • 34 + 19 • 1^5 - 21 • 25 - 6 • 16 + 7, 

/ (13^) = -552 • 1^2^ + 748 • 12* 

-406 • 1*23 + 2622 • 1^223 + 500 • 1^3^ - 2180 • 123^ 

-1226 • 2^3 + 280 • 3^ - 26 • 1^4 + 720 • 1^24 - 965 • 1^34 - 1620 • 12^4 + 280 • 14^ 

+1260 • 234 + 67 • 1*5 - 576 • 1^25 + 375 • 2^5 + 445 • 135 - 126 • 45 - 53 • 1^6 

+221 • 126 - 84 • 36 + 28 • 1^7 - 36 • 27 - 8 • 18 + 9, 

/(123) = 10 • 12^ + 5 ■ 1^3 - 10 • 23 - 4 • 14 + 5, 

/ (1232) ^ _72 . i223 _ 48 • 12^ + 48 • 13^ 

+90 • 2^3 - 8 • 1^4 + 68 • 124 - 35 • 34 + 18 • 1^5 

-20 • 25 - 5 ■ 16 + 7, 

/(133) = -6- 23-2- 14 + 5. 
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We now give /s, but only when = 0, as the general case is too long to include here: 

r = 5 : /(7) = Hq is given by (|R5|) . 

/(45) ^ -19 • 1^2 + 189 • l'^2^ - 411 • 1^2^ + 105 • 2'^ 

+51 • 1^3 - 542 • 1^23 + 276 • 1^3^ + 837 • 12^3 - 280 • 23^ 

-69 • l'^4 + 417 • 1^24 - 210 • 2^4 - 279 ■ 134 + 35-42 

-168 • 125 + 56 • 1^5 + 56 • 35 - 28 • 1^6 + 28 • 26 + 8 ■ 17 - 8, 

/(25) = -9 • 1^2 + 42 • 1^2^ - 15 • 2^ + 16 • 1^3 

-60 • 123 + 10 ■ 3^ - 14 • 1^4 + 15 ■ 24 + 2 ■ 15 - 6, 

/(16) = -5 • 1''2 + 30 • 1^2^ - 15 • 2^ + 10 • 1^3 

-50 • 123 + 10 • 3^ - 10 • 1^4 + 15-24 + 15-6, 

/ (14^) ^ 75 • 1^2^ - 288 • 1^2^ + 105 ■ 2^ 

+15 • 1^3 - 324 • 1^23 - 33 • l'^4 + 723 • 12^3 + 204 • 1^3^ 

-279 • 23^ + 309 • 1^24 - 210 • 2^4 - 243 • 134 + 35-42 

+35 -1^5 - 147 -125 + 56 -35 

-21 -126 + 28 -26 + 7- 17 - 8, 

/(124) = 21 • 1^2^ - 15 • 2^ + 7 • 1^3 - 48 • 123 

+10 ■ 32 - 9 • 1^4 + 15-24+15-6, 

/(122) = 3-22 + 3 - 13-4, 

/ (1^5) = 12 • 1^2^ - 12 • 2^ + 4 • 1^3 - 36 • 123 + 10 • 32 - 6 ■ 1^4 + 14-24-6, 
/ (1^4) = -6 ■ 2^ - 18 • 123 + 9 • 32 - 3 ■ 1^4 + 12 • 24 - 15 - 6, 
/(l'2)^-4. 

We skip /g. The extra g{iT) needed for gr, r = 2, 3, 4 in terms of a as follows: 
r = 2 : 

g (32) = -4 • 1^2 + 10 • 12^ + 8 • 1^3 - 8 • 23 - 5 • 14 + 5. 
r = 3 : 

5(34) = -6 • l'*2 + 30 • 1^2^ - 12 • 2^ + 15 • 1^3 

-51 • 123 + 9 • 32 + 14 • 24 - 14 ■ 1^4 + 6 • 15 - 6, 

g (3^) ^ 48 • l'^22 - 212 • 1^2^ + 68-2'' 

-284 • 1^23 + 8 ■ 1^3 + 594 • 12^3 + 226 • 1^3^ - 226 • 23^ 

-28 • l'*4 + 306 • 1^24 - 168 • 2^4 - 265 • 134 + 35-42 

-144 -125 + 38 -1^5 + 53 -35 

-25 -126 + 25 -26 + 8 -17- 8, 

5(23) = -2 - l22 + 2 - 22 + 3 - 13 - 4. 
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r = 4 : 

5(35) = -8 • 1^2 + 68 • 1^2^ - 84 • 12^ 

+24 • 1^3 - 176 • 1^23 + 66 • 13^ + 92 • 2^3 - 30 • 1^4 

+98 • 124 - 34 • 34 + 20 • 1^5 - 20 • 25 - 7 • 16 + 7, 

g (3^4) = 84 ■ 1^2^ - 618 • 1^2^ + 642 • 12"* 

-600 • l'*23 + 2652 • l22^3 + 12 • 1^3 

-2136 ■ 123^ - 1002 • 2^3 + 624 • 1^3^ + 252 • 3^ 

-48 • 1^4 + 863 • 1^24 - 1565 ■ 12^4 

-1126 • 1^34 + 1141 • 234 + 310 • 14^ + 79 ■ l'*5 - 629 • 1^25 + 334 • 2^5 + 483 • 135 - 125 • 45 
-69 • 1^6 + 231 • 126 - 81 ■ 36 + 34 • 1^7 - 34 • 27 - 9 • 18 + 9, 
g (3*) = -1008 • 1^2^ + 6752 • 1^2'^ - 6044 • 12^ 

-416 • 1^23 + 12224 • l'*223 - 41344 • 1^233 + 12124 • 2'*3 + 1264 • 1^3^ 
-27728 • 1^232 + 52552 • 12^3^ + 12896 • 1^3^ - 12896 ■ 23^ 

-16 • 1^4 + 1840 ■ 1^24 - 19320 • 1^224 - 15395 • 124^ - 5864 • l''34 + 53728 • 1^234 
-28524 • 2^34 - 21740 ■ 13^4 + 25276 • 12^4 + 4200 • 1^4^ + 5635 • 34^ 
+80 • 1^5 - 3320 • 1''25 + 14988 • 1^225 - 5480 • 2^5 

-24152 • 1235 + 6760 ■ 1^35 + 4348 • 3^5 - 6198 ■ 1^45 + 6348 • 245 + 1353 • 15^ 
-168 • 1^6 + 3136 • 1^26 - 5704 • 12^6 - 4136 ■ 1^36 
+4136 • 236 + 2280 • 146 - 456 ■ 56 

+192 • 1''7 - 1640 • 1^27 + 852 • 2^7 + 1288 • 137 - 330 • 47 
-129 ■ 1^8 + 451 • 128 - 161 • 38 
+51 • 1^9 - 51 • 29 - 11 • 1(10) + (11), 
g (23^) = 20 ■ 1^2^ - 50 • 12^ + 4 • l''3 
-96 • 1^23 + 74 -223 + 54 -13^ 

-12 • 1^4 + 73 • 124 - 33 • 34 + 13 • 1^5 - 18 • 25 - 6 • 16 + 7. 



Appendix E: Expansions of fr, Qr for the gamma when /a = 

Here, we give the coefficients /(vr), ^'(vr) needed in (I1.34D for /r, gr-, r < 4: oi h when ^3 = 
for X = G a, gamma variable with mean m in terms of a = m — 1, y = — l/y, as in (j2.2p . 

One can show 

/ (1\ fc + 2) = {-l)\k + l)c7'^+^+i ^ H,+ir+^+i [j + 2),_i, ^ > 1, (E.l) 
/ {l\ 2, + 1) = (-1)- ^ O'^ {(2fc + l)(a - j) + f) {j + 2),. (E.2) 

The other /(vr) are more simply given by the formulas for them in Appendix C, but we 
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give them here for the record: 

r = 2 : for fc = 4, 2, /(fc) ^ of (lO)) . 

r = 3 : /(tt) of ((EJ|) for tt = 5, 12, 14. 

r = 4 : /(tt) of jEl]), jE!]) for tt = 6, 15, 24, 2^ 1^4, 1^2, 25, 
/ (4^) /3o-^ = -5af2 - a(25a - Tl)-f - 10[a]2(5a - l)t 
-2[q;]2 {iho? ~ 89a + 84) - Ha {^ho? ~ 109a + 140) 
-Ha (Sa^ - 39a2 + 114a - 120) /. 

r = 5 : /(tt) of ([EI]), for tt = 7, 16, 25, 124, 12^, 1^5, l34, 1^2, 

/(45)/cr^ = -19a/ _ 6a(19a - 17)/ - 3[a]2(95a - 138)/ 

-4[a]2 (95a2 - 349a + 336) / - 3[a]3 (95a2 - 433a + 560) / 

-6[a]3 (I9a3 - 154a2 + 455a - 480) / - \a\^ (l9a^ - 177a^ + 638a - 840) /, 

/ (14^) /3cr* = 10a/ + 3a(25a - 22)/ + 40[a]2(5a - 7)/ 

+10[a]2 (25a2 - 89a + 84) / + 6[a]3 (25a2 - 109a + 140) / 

+7[a]3 (5a^ - 39a2 + 114a - 120) /. 
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r = 6 : /(tt) of dEI]), dElJ for tt = 8, 17, 26, 1^6, 1^5, 1*4, 1*2, 125, 1^24, 1^2^, 
/ (5^) /2cr^ = -I2af - 4a(21a - 19)y^ - 12[a]2(21a - 31)/ 
-12[a]2 (35^2 - 130a + 126) / 

-20[a]3 (21^2 - 97a + 126) / - 12[a]3 {2lo? - Y12o? + 511a - 540) / 

-12[a]4 (7a^ - 66^^ + 239a - 315) - 4[a]4 (3a'^ - 43a3 + 258a2 - 743a + 840) , 

/(46)/cr^ = -23a/ - a(161a - 146)/ - 3[a]2(161a - 242)/ 

-5[a]2 (I61a2 - 610a + 600) / 

-5[a]3 (laia^ - 767a + 1008) / 

-3[a]3 (laia^ - 1357a2 + 4082a - 4320) / 

-[a]4 (I61a3 - 1575a2 + 5734a - 7560)/ 

-[a]5 (23a3 - 249a2 + 1066a - 1680) /, 

/ (4^) /9o-" = 14a/ + 3a(57a - 50)/ 

+a (805a2 - 1860a + 1064) / 

+ [a]2 (2023a2 - 6620a + 5880) / 

+3[a]2 (10150^ - 5840a2 + 12056a - 8720) / 

+ [a]2 (2849a'' - 24931a3 + 87918a2 - 145044a + 92400) / 

+ [a]3 (1631a'' - 16737a^ + 71986a2 - 149400a + 123200) / 

+ [a]3 (525a5 - 7500a* + 47457a^ - 162474a2 + 294472a - 221760) /° 

+ [a]3 (73a^ - 13740^ + 11877a* - 591480^ + 1753720^ - 288080a + 201600) /\ 

/ (24^) /cr^ = 62a/ + 81a(7a - 6)/ 

+3a (635a2 - 1450a + 824) / 

+22[a]2 (I45a2 - 473a + 420) / 

+36[a]2 (mo? - mZo? + 961a - 700) / 

+3[a]3 (449a3 - 3079a2 + 8006a - 7560) / 

+ [a]3 (257a* - 26910^ + 117460^ - 24552a + 20160) /, 

/ (2^4) /3cr^ = 10a/ + a(61a - 50)/ 

+a(l23a2- 272a + 152)/ 

+ [a]2 (I03a2 - 326a + 280) / 

+ [a]3 (31a2 - 114a + 120) /, 

/ (2^) /cr^ = 14a/ + a (43a - 30)/ 

+a (29a2-50a + 24)/, 

/(145)/2o-^ = 19a/ + 9a(19a - 17)/ + 6[a]2(95a - 138)/ 

+10[a]2 (95a^ - 349a + 336) / 

+9[a]3 (95a2 - 433a + 560) / 

+21[a]3 (I9a3 - 154a2 + 455a - 480) / 

+4[a]4 (I9a^ - 177a^ + 638a - 840) /, 

/ (1^4^) /6cr^ = -15a/ - 6a(25a - 22)/ - 100[a]2(5a - 7)/ 
-30[a]2 (25a2 - 89a + 84) / 
-21[a]3 (25a2 - 109a + 140)/ 
-28[a]3 (5a^ - mc? + 114a - 120) /. 
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When ^3 = 0, g{Tr) needed for gj- in (II .341) are as follows: 



g(2, fc + 1) = ^ H,+i {ka - (fc - (E.3) 

r = 2: for A: = 4, 2, ^(/s) = Hk-i of (lO]) . 
r = 3 : 5(5) = of (1231) . 

r = 4 : 5(6) = iJg of 1^ and 5(24),5 (22) of (Ie3]) . 

g (4^) /3cr^ = -3a/ - 3a(5a ~ 6)j/3 - 6H2(5a - 11)/ 

-6[a]2 (Sa^ - 25a + 28) / - [a]^ {l5a^ - 99a + 140) / 

-[a]3(a - 6) (3a2 - 15a + 20) f. 

r ^ 5 : 5(7) = He of ([231), 5(25) of (jEH, 

5(45)/12cr* = -a/ - a(6a - 7)/ - [a]2(15a - 32)/ 

-2[a]2 (lOa^ - 49a + 55) / - H3 (l5a2 - 97a + 140) / 

-[a]3 (6a3 - 65a2 + 219a - 240) / - [a]4(a - 7) (a^ - 6a + 10) /. 



r = 6 : 5(8) = Hr of 5(26) of (Ie3| . 

5 (5^) /8cr^ = -2a/ - 2a(7a - 8)/ - 3[a]2(14a - 29)/ - [a]2 {70a^ - 335a + 372) / 

-10[a]3 (7a2 - 44a + 63) / - 6H3 {7a^ - l^o? + 247a - 270) / 

-[a]4 (I4a^ - 177a2 + 703a - 945) / - [a]4(a - 8)(2a - 7) {o? -7a + 15) /, 

5(46)/15cr^ = -a/ - a(7a - 8)/ - [a]2(21a - 44)/ 

-H2 (35a2 - 170a + 192) / - [a]3 (35a2 - 225a + 332) / 

-[a]3 (21a3 - 227a2 + 776a - 864) / - [a]4 {7a^ - 91a^ + 372a - 504) / 

-[a]5(a-8) (a2-7a+14)/, 

5 (43) /9cr" = 6a/ + 12a(5a - 6)/ + 3a (liga^ - 380a + 264) / 

+3[a]2 (301a2 - 1480a + I68O) / + 3[a]2 (455a3 - 3800a2 + 9864a - 8I6O) / 

+3[a]2 (427a* - 5333a3 + 234540^ - 44092a + 30240) / 

+ [a]3 (735a4 - 11265a3 + 60786a2 - 141816a + 123200) / 

+ [a]3 (237a^ - 4944a* + 39033a3 - 1505580^ + 288712a - 221760) /° 

+ [a]3(a - 10) (33a5 - 564a* + 3957a3 - 142980^ + 26552a - 20160) /\ 

5 (24^) /3ct^ = 6a/ + 3a(19a - 24)/ + 3a (65a2 - 222a + I6O) / 

+6[a]2 (55a2 - 297a + 360) / + 12[a]2 (25a3 - 2300^ + 642a - 560) / 

+ [a]3 (Mla^ - 1665a2 + 5948a - 6720) / 

+ [a]3(a - 8)(3a - 10) {9a^ - 53a + 84) /, 

5 (2^4) /3ct^ = 2a/ + a(13a - 18)/ + a {27a^ ~ 104a + 80) / 

+ [a]3(23a - 100)/ + [a]3(a - 6)(7a - 20)/, 

g (2^) /o-s = 2a/ + a(7a - 12)/ + a(a - 4)(5a - 6)/. 
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Appendix F: Comparison with Cornish and Fisher 

Here, we specialize to the case X = A^, so that = Her, the rth Hermite polynomial of 
(fL36]) . We give /(vr), g{7r) of (fLMD needed for fr, Qr of (fTTOD and (fTTT]) . We also give 
special formulas for some of them that do not hold in the non-normal case. Recall that for 
z > 2, E A^^-?' = 1 • 3 • • • (2j - 1). Using the notation (fOHjl . we have 

/(l^A;)=0, f{V,k + l) = {-iy[k],Hk-,, f{2^) = {-iy-'H,EN'^. 

As above, special relations are put in the line starting 'r ='. /(vr) needed for fr are as 
follows: 

r = l: f{3)^H2, /(I) = 1- 

r = 2 : /(4) = iJg, /(2) = H,, /(13) = -2ffi, 

/ (32) = -2 (4a;3 _ 7a;) . 

r = 3 : /(5) = H^, /(14) = -3il2, /(12) = -1, / (l^S) = 2, 
/(34) = -llx" + 42a;-'^ - 15, / (3^) = 2 (69a;'' - 187a;^ + 52) , 
/(23) ^ -5a;2 + 3, / (13^) = 2 (l2a;2 - 7) . 

r = 4 : /(6) = H^., f (2^) = -3i/i, /(15) = f (1^4) = 6i7i, 

/(l22) =/(l33) -0, 

/ (42) = -3 (5a;^ - 32x^ + 35a;) , /(35) = -2 (7a;^ - 48a;^ + 51a;) , 
/ (3^4) = 2 (llla;^ - 547a;3 + 456a;) , / (3'') = -4 (948a;-^ - 3628a;3 + 2473a;) , 
/(24) = - (7a;3 - 15a;) , / (23^) ^ 2 {36x^ - 49a;) , /(134) = 4 (lla;'"^ - 21a;) , 
/ (13^) = -4 (I38a;=' - 187a;) , /(123) = lOHi, f (1^3^) = -A8H1. 



r = 5 : /(7) = He, /(16) = -5^4, / (I's) = I2H2, f (iH) = -6, 
/(l32) -/(143) =0, 

/(45) = -19x^ + 189a;^ - 411a;2 + 105, /(36) = -17x^ + 185a;^ - 405a;2 + 105, 

/ (342) = 3473;*^ - 2643a;'' + 45213;^ - 945, / (3^5) = 2 (l62a;S - 1309a;4 + 2232a;2 - 471) , 

/ (3^4) = -2 (3354a;^ - 20831a;'' + 29148a;2 - 5174) , 

/ (3^) = 4 (36240x^ - 184146x* + 217921a;2 - 33523) , 

/(25) = -3 {3x^ - 14a;2 + 5) , /(234) = 1213;" - 378x'^ + 105, 

/ (23^) -2 (8973;" - 2057a;2 + 468) , / (2^3) 5 {7x^ - 3) , 

/ (142) = 3 (25a;4 - 96a;2 + 35) , /(135) = 2 {35x'^ - lUx^ + 5l) , 

/ (13^4) = -6 (l85a;^ - 547a;2 + 152) , / (13") = 4 (4740^" - 10884a;2 + 2473) , 

/(124) = 3 {7x^ - 5) , / (123^) = -2 (l08.x2 - 49) , / (12^) ^ 3, 

/ (1^34) = -12 (lla;2 - 7) , / (1^3^) = 4 (414.a;2 - 187) , / (1^23) = -10, / (1^3^) = 48. 
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/e has forty six terms, so is only given here for the case I3 = 0: 

r = 6 : /(8) = i/7, /(17) = -6^5, / (l'6) - 20H3, f (l^S) - f (2^) = 15ffi, 

/ (l^S) = / (1^4) = / (l^S) = / (iH) = f (iH^) = / (1323) = 0, 

/ (52) = -24 {x^ ~ Ux^ + 51x^ - 39x) , /(46) = ~23x^ + 333x^ ~ 1215x^ + 945x, 

/ (4^) = 3 (l77x^ - 1899x^ + 54512;^ - 3465a;) , /(26) = -llx^ + 90^^ - 105a;, 

/ (24^) = 3 (65a;^ - 352a;3 + 315a;) , / (2^4) = 21 {3x^ - 5x) , 

/(145) = 6 (I9a;5 - 126a;3 + 137x) , /(125) = 12 {3x^ - 7x) , 

/ (1^42) = -12 (25a;3 - 48x) , / (1^24) = -42Hi. 

Note the special relations 

/ (1232) = -A8H1, f (1^24) = -42Hi, /(123) = lOHi. 

For A: > 2, i > 0, gi2\k) = {-Ifi^kzHk-u where i^ko = 1, i^ki = {k-l){k + l) ■ ■ ■ {k + 2i-3) 
for i > 1. That is, 

g {2\r + 2-2i) = {-If < r >i F,+i_2*, < i < r/2, 

where 

< r >o= 1, 

< r >i= [r - l)(r - 3) • • • (r - 2i + 1) = 2'T(s)/r(s - i), s = (r + l)/2, i > 1. 

Again, cases of this relation are put in the line next to 'r ='. The coefficients (/(vr) of (ll.34p 
needed for are as follows: 

r = l: 9(3) =H2, .9(1) = 1. 
r-2: .9(4)=i/3, 9(2) = H,, 
g (3^) = -2 (2a;3 - 5a;) . 
r = 3 ; .9(5) - i/4, 5(23) - -2i/2, 

g(34) = -6 (x-* - 5x2 ^ 2) , g (33) = 4 (12a;'' - 53^^ + 17) . 

r = 4 : 5(6) = i/5, .9(24) = -3i?3, .9 (2') = -H,, 

g (42) = -3 (3a;5 - 24a;3 + 29x) , g(35) -4 (2a;'^ - 17a;3 + 21a;) , 

g (3^4) = 6 (I4x^ - 103a;3 + 107a;) , 9 (S'^) = -4 (252a;^ - leSSa;^ + 1511a;) , 

g (232) = -55 (32) . 

r = 5 : 5(7) = H^, .9(25) = -4if4, .9 (2^3) - Sila, 

5,(45) = -12 (a;*^ - 12x'* + 29a;2 - S) , ^(36) = -10 (x^ - 13a;'' + 33^^ - 9) , 

g (34^) = 12 {I2x^ - 129a;'' + 271a;2 - 64) , .9 (3^5) = 8 (l6x^ - 181x'' + 393^^ - 90) , 

g (334) = -24 (80x^ - 803a;'' + 1513^^ - 304) , 

g (3^) = 32 (960x6 _ 89373,4 ^ i^qq2x^ „ 2651) , 

5(234) = -65(34), 9 (233) = -85 (33) . 
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r = 6 : 5(8) = F7, 5(26) = -bH^, g (2^4) = I5H3, g {2') = BH^, 

g (5^) -8 {2x^ - 33x-^ + U2x^ - 108x) , 

.g(46) = -15 (x^ - 17x'^ + 69x'^ - 57x) , 

g (4'^) = 27 (9x^ - 131a;^ + 451a;^ - 321a;) , 

g{S7) = 6 (2a;'' - 37a;^ + 160a;^ - 135a;) , 

5(345) = 12 (18a;'' - 273a;^ + 974a;^ - 695a;) , 

g (3^6) = 10 (I8x'' - 2933;^^ + llOOa;^ - 795a;) , 

g (3^4^) = -6 (594a;'' - 8193a;^ + 26006x^ - 16367a;) , 

g (3^5) = -8 (396x'' - 5708x^ + 18755a;^ - 11811a;) , 

g (3^4) = 12 (5148a;'' - 67004.x''^ + 195259x^ - 109553x) , 

g (36) = -8 (I54440x'' - 1887684.T-^ + 5033714x^ - 2542637x) , 

g (242) = -7g (42) , 5(235) = -7g(35), g (23^4) = -9g (3^4) , 

g(234)=-llg(34), 5(2^3^)= 355(32). 

These results agree with those given on page 317 of Cornish and Fisher (1937) and page 
214 of Fisher and Cornish (1960), except for (i) a typo on page 316: ab was written for ad, 
making /(14) above into a second /(12); and (ii) they gave /(15) = — 4(x^ — x) instead 
of —AH3. The factor —6 in 5(34) was omitted by Hill and Davis (1968). x in /(3^) was 
omitted in Withers (1984). Fisher and Cornish did not give h,., /,• for r = 5,6. Hill and 
Davis reference a Stanford Technical Report which may give er{x, L) for e = h, f , g up to 
r = 11 when X = N. 

Note how 5(27r) is a multiple of ^(Tr). For example, if tt = (3*^4*4 • • • ), then g{2'K) = 
(1 - k|)5(7r), where |7r| = Eit=3^^fe- 
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Table 1.1 Number of terms needed by er{x) for different choices of X. 
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Table 2.1 Two approximations for the 95th quantile of 2 ^ lnF24,60- 



Order of 
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totals 


errors 
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Table 2.2 Two approximations for the 95 quantile of 2 ^ InFs 5. 



Order of 


Successive 


Successive 


Successive 


magnitude 


terms 


totals 


errors 





.0 






1 


-.0 




-.00 


2 


.0 




.00 


3 


-.0 




-.00 


4 


.0 




.00 


5 


-.0 




-.00 


6 


-.0 




.00 



43 



